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ABSTRACT 

This  project  was  undertaken  to  study  earthquakes  that  occur  in  volcanic 
regions  and  which  result  from  volcanic  processes .  The  purpose  was  to  deter¬ 
mine  whether  or  not  these  volcanic  quakes  are  the  same  as  ordinary  eartht- 
quakes,  and  if  not,  just  how  they  differ  and  what  seismic  parameters  can  be 
used  to  identify  them. 

Three  areas  were  selected  for  investigation:  Italy,  Hawaii  and  Japan. 
Volcanoes  in  Italy  and  Japan  are  explosive  while  those  in  Hawaii  are  not. 

The  Italian  volcanoes  have  been  under  investigation  for  many  years  and  under 
observation  for  centuries.  Hawaiian  Volcanology  and  Japanese  Volcanology 
are  more  recent  but  both  have  contributed  greatly  to  the  application  of 
modem  geophysical  techniques  to  Volcanology. 

Seismic  data  from  all  three  areas  were  analyzed.  Crustal  models  and 
seismic  wave  velocities,  where  available,  were  used  as  a  basis  for  locating 
epicenters.  Graphical  methods  such  as  the  A  P-H  Arc -Distance  Method  and  the 
Hyperbolic  Method  were  used  to  locate  epicenters.  A  computer  program  was 
set  up  for  the  same  purpose.  Several  techniques  were  developed  to  augment 
the  use  of  data  in  the  computer  program.  One  of  these  was  the  technique  of 
supplying  a  fictitious  arrival  time  at  a  seismic  station  and  varying  this 
arrival  time  value  until  all  the  input  data  converged  to  a  solution.  Results 
of  the  graphical  solutions  and  those  from  the  computer  program  were  in  excel¬ 


lent  agreement. 


Epicentral  locations  gave  distances,  azimuths  and  foc&  depths  vnich 
were  important  in  analyzing  various  seismic  parameters.  Hawaiian  volcanic 
quakes  show  distinct  tectonic  characteristics.  Quadrantai  patterns  of  first 
motion  can  be  delineated  and  related  to  tectonic  features.  Japanese  quakes 
show  a  mixed  set  of  characteristics.  Some  patterns  are  quadrantai,  indicat¬ 
ing  a  tectonic  relationship,  while  others  show  push  or  compres sional  motion 
at  all  stations  indicating  an  explosive  mechanism.  Vesuvius  shows  b«.'„h  com¬ 
press  ional  and  dilatat ional  first  motion. 

Amplitude  distance  relationships  fail  to  show  statistical  significance 
at  Vesuvius  because  the  observations  are  all  from  a  single  station.  In 
Hawaii  there  is  a  significant  relationship.  Magnitude  curves  based  on  log  A 
versus  log  4  show  that  the  earthquakes  are  of  low  energy  with  magnitudes 
from  plus  one,  down  to  negative  values  and  quakes  are  usually  not  recorded 
at  any  great  distance.  Attenuation  is  high  in  the  caldera  region. 

Volcanic  quakes  occur  in  swarms  and  an  analysis  of  the  frequency  of  oc¬ 
currence  was  made.  Tectonic  and  explosion  quakes  have  characteristic  values 
for  m  in  the  equation  NAm  =  c.  This  parameter  or  coefficient  may  depend  or. 
a  number  of  factors  -  such  as  focal  mechanism,  focal  depth,  and  stress  ap¬ 
plication.  Stress  application  certainly  is  part  of  the  focal  mechanism  but, 
in  the  sense  of  an  externally  applied  surface  stress,  it  is  vastly  different 
from  any  of  the  usual  earthquake  mechanisms.  It  seems  possible  that  such  a 
stress  application  might  have  a  characteristic  m  value. 


ANALYSIS  OF  THE  RELATIONSHIP  BETWEEN 
ASSOCIATED  VOLCANIC  AND  SEISMIC  EVENTS 


1 

INTRODUCTION 

Perhaps  one  of  the  oldest  associations  of  natural  phenomena  is  that  of 
volcanoes  and  earthquakes.  In  the  many  regions  of  the  earth  In  which  these 
phenomena  occur,  man  has  puzzled  as  to  their  relationship.  In  more  modern 
times  geophysicists  have  come  to  a  recognition  of  tectonic  earthquakes  as 
distinct  from  volcanic  earthquakes  but  even  here  the  distinction  becomes 
vague  in  many  cases. 

Volcanic  earthquakes  could  be  loosely  described  as  earthquakes  assoc¬ 
iated  with  a  volcanic  eruption,  leaving  -unanswered  many  questions  as  to  the 
seismic  events  in  both  the  pre -eruptive  and  post -eruptive  periods.  One  set 
of  notable  characteristics  of  volcanic  earthquakes  is  that  they  are  very 
low  in  intensity  and  occur  in  swarms. 

The  nature  of  these  seismic  events  can  be  better  perceived  when  the 
fundamental  volcanic  process  is  known  and  understood.  The  cause  and  me¬ 
chanics  of  a  volcanic  eruption  represent  a  very  complex  subject  which  the 
author  prefers  to  leave  to  those  who  are  better  qualified.  However,  all 
eruptions  involve  certain  basic  processes  which  can  b"  delineated  here  and 
which  will  throw  some  light  cn  the  resulting  seismic  events.  First,  there 


is  t,he  pre-eruptive  stage  during  which  there  is  the  movement  of  magmatic  ma¬ 
terial  from  depth  toward  the  surface  of  the  earth.  Associated  with  the  ver¬ 
tical  rise  of  the  magma  there  may  be  significant  horizontal  movements  as 
well.  Second,  tt-r''  is  the  actual  eruptive  stage  in  which  volcanic  material 
is  discharged  at  the  surface.  Tne  eruptive  stage  may  vary  from  a  quiescent 
outpouring  of  fluid  lava  at  the  one  extreme  to  the  violently  explosive  erup¬ 
tion  at  the  other  extreme.  Volcanologists  have  developed  an  impressive 
classification  of  eruptive  types  and  resulting  eruptive  materials.  Finally, 
there  is  the  post-eruptive  stage  during  which  the  processes  that  have  been 
active  in  producing  the  eruption  are  seeking  an  equilibrium  condition. 

During  these  various  stages,  seismic  activity  manifests  itself  as  an 
increasing  swarm  of  very  small  earthquakes.  These  are  probably  associated 
with  the  transport  of  magma  either  vertically  or  horizontally,  resulting  in 
the  splitting  and  rupture  of  the  crustal  rock  layers  and  with  the  general 
distention  of  the  area.  Closely  associated  with  the  actual  eruptive  stage 
is  the  development  of  volcanic  tremor  which  is  a  continuous  vibration  of  the 
crust. 

2 

AREAS  INVESTIGATED 

The  present  project  has  as  its  immediate  purpose  the  investigation  of 
seismic  signals  originating  in  volcanic  areas  to  delineate  if  possible 
unique  identifying  charar  .eristics  which  will  distinguish  them  from  seismic 
signals  associated  with  other  types  of  earthquakes  whether  natural  or  arti¬ 
ficial.  Originally  two  volcanic  areas,  Italy  and  Hawaii,  were  selected  for 
investigation,  and  then  a  third  area,  Japan,  was  added. 


3 

MOUNT  VESUVIUS 


The  Italian  investigation  began  with  a  visit  to  Professor  Pietro  Caloi, 
Director  of  the  National  Institute  of  Geophysics.  Here  data  and  information 
on  the  local  crustal  structure  and  seismicity  were  gathered  as  background 
information  for  the  investigations  of  Mount  Vesuvius  and  Mount  Etna  which 
were  to  follow. 

3.1  THE  VESUVIAN  OBSERVATORY 

At  Mount  Vesuvius,  Professor  Guiseppe  Imbo,  Director  of  the  Vesuvian 
Observatory  and  also  Director  of  the  Institute  of  Physics  c  '  the  Earth, 
University  of  Naples,  put  the  facilities  of  the  Observa-  cry  at  my  disposal. 
Due  to  the  lack  of  instrumental  constants  only  data  after  1944  were  useful 
for  the  investigation.  The  instruments  consisted  of  two  horizontal  com¬ 
ponent  200  kg  W:\echert  Seismographs  and  a  vertical  component  80  kg  Wiechert 
Seismograph. 

Free  period  and  static  magnification  for  the  instruments  are  as  follows: 
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4.1 

158 
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Magnification  curves  were  available  for  each  instrument.  The  approximate 
distance  from  the  Vesuvian  Observatory  to  the  center  of  the  crater  of  Mount 
Vesuvius  is  2.5  kilometers.  The  elevation  of  the  Observatory  is  609  meters, 
that  of  the  crater  floor  is  951  meters,  and  the  rim  of  the  crater  is  some 
200  meters  higher.  This  presents  the  interesting  situation  that  seismic  foci 
could  be  either  above  or  below  the  Observatory  level. 
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3.2  VESUVIAL  EARTHQUAKES 


In  the  following  table,  time  is  Central  E  opean  Mean  Time,  Amplitude  A 
is  in  millimeters  as  read  from  the  record,  period  T  is  in  seconds,  and  direc 
tion  of  motion  is  N-north,  S-south,  E-east,  W-west,  c-compression,  and  d- 
dilatation. 


TABLE  I  VESUVIAN  QUAKES 


No. 

Date 

Phase 

Time 

Direction 

Amplitude 

Period 

1 

11-3-61 

iP 

18-28-25.0 

S 

0.2 

iP 

25.0 

- 

- 

- 

iP 

25.6 

d 

0.1 

- 

eS 

25.9 

N 

0.7 

- 

eS 

26.0 

E 

0.2 

0.2 

2 

8-4-61 

iP 

01-50-27.7 

S 

0.4 

0.4 

iP 

28.1 

E 

0.3 

0.4 

3 

7-27-61 

eP 

13-01-36.3 

N 

0.2 

0.4 

iP 

37.1 

W 

0.4 

0.6 

4 

7-6-61 

iP 

11-22-23.6 

s 

0.4 

0.4 

eP 

23.6 

E 

0.4 

0,6 

eS 

25.2 

E 

0.7 

0.8 

iS 

25.4 

S 

1.0 

0.5 

5 

6-27-61 

iP 

05-02-08.3 

S 

0.2 

0.3 

eP 

08.3 

E 

0-3 

0.6 

6 

5-19-61 

iP 

23-05-42.5 

d 

— 

- 

iP 

42.6 

N 

0.2 

0.6 

e? 

43.0 

W 

0.4 

0.6 

7 

5-11-61 

iP 

23-53-28.,. 

s 

0.3 

0.5 

iP 

28.5 

E 

0.4 

0.5 

IS 

29.0 

N 

0.8 

0.7 

iS 

30.0 

E 

0,8 

0.5 

8 

5-10-61 

iP 

17-16-38.8 

E 

0.3 

0.1 

e 

41.7 

W 

0.7 

0.4 

i 

42.6 

s 

1.0 

0.6 

4 


No. 

Date 

Phase 

Time 

9 

9-29-60 

IP 

10-32-11.8 

IP 

11-9 

iP 

12.0 

iS 

12.6 

IS 

12.8 

10 

8-30-60 

iP 

07-3^-33-7 

iP 

34. 4 

iS 

3  b.l 

iS 

34.7 

11 

7-1U-60 

iP 

04-00-35-3 

iP 

35-1 

iP 

35-4 

12 

7-  7-60 

eP 

01-16-32.1 

eP 

32.1 

13 

7-  7-60 

IP 

OI-I0-33.5 

iP 

33-9 

iP 

34.0 

IS 

37.8 

i 

39-1 

14 

12-  2-56 

eP 

02-54-55-2 

iP 

56.4 

IP 

56.4' 

iP 

56.4 

iS 

58.0 

IS 

58.3 

iS 

58.8 

i 

55-00.2 

15 

12-  2-56 

iP 

02-42-55-2 

eP 

55-4 

iP 

55-5 

iS 

58.3 

iS 

58.6 

i 

43-01.1 

16 

11-30-56 

iP 

22-33-41.7 

iP 

42-4 

iP 

42.4 

iS 

44.5 

iS 

44,6 

iS 

45.2 

i 

46.4 

Direction 

Amplitude 

Period 

N 

0.2 

0.2 

d 

0.2 

0.2 

W 

0,4 

0,2 

S 

1.0 

0.4 

E 

o.;n 

0.2 

W 

0.3 

0.4 

d 

0.3 

0.2 

N 

0.8 

0-3 

E 

0.8 

0.6 

S 

0.1 

- 

W 

0  2 

0.2 

d 

0.1 

- 

S 

0.1 

0.8 

- 

0.1 

0.2 

N 

0.2 

0.6 

W 

0.2 

0.5 

a 

0.2 

0.4 

s 

0.5 

0.7 

E 

0.5 

0.5 

d 

0.3 

0.8 

c 

0.5 

0.8 

N 

0.3 

0.8 

W 

0.4 

n.8 

W 

2,2 

0.8 

S 

2.6 

0.9 

c 

0.7 

0.8 

E 

2.7 

0.9 

W 

0.3 

0.6 

c 

0.3 

0.6 

N 

0.3 

0.4 

d 

0.4 

0.6 

N 

0.6 

- 

•"i 

O 

0.8 

1.0 

E 

C.7 

0.8 

d 

0.5 

0.8 

N 

1.6 

0.8 

c 

0.7 

0.4 

E 

2.8 

0.6 

N 

3-3 

0.6 

S 

5-4 

0.6 

5 


No. 

Tate 

Phase 

Tim- 

Direction 

Amplitude 

Period 

IT 

11-30-56 

IP 

23-33-09.8 

S 

0.5 

0.8 

/X  s' 

•?  TO 

X X 

09.9 

W 

0.5 

0.6 

IP 

10.4 

d 

0.2 

0.5 

IS 

13.8 

d 

0.5 

0.6 

iS 

13.8 

N 

0.5 

0.6 

iS 

14.2 

W 

0.8 

0.6 

18 

9-24-56 

iP 

22-40-45-7 

d 

0.3 

0.4 

iP 

45.9 

S 

0.8 

0.6 

eP 

45.9 

E 

0.6 

0.6 

iS 

46-9 

S 

0.2 

1.6 

iS 

47-1 

w 

0-7 

0.8 

eS 

47.6 

d 

0.8 

0.4 

19 

7-26-56 

eP 

03-03-42.5 

d 

0.2 

- 

iP 

42.6 

S 

0.4 

— 

iP 

42.6 

Hi 

0.4 

— 

20 

3-22-56 

IP 

18-06-39.6 

W 

0.3 

0.2 

iP 

39.8 

d 

0  •  i. 

0.4 

iP 

39.8 

S 

0.4 

0.4 

iS 

40.2 

W 

0.3 

- 

iS 

40.7 

d 

0.7 

0.4 

iS 

40.7 

N 

0.3 

0.4 

i 

41.3 

S 

0,8 

1.6 

21 

2-2-56 

eP 

16-53-23.2 

W 

0.4 

- 

IP 

23.5 

c 

0.2 

0.2 

IP 

25.0 

c 

0.2 

0.2 

eP 

25.7 

N 

0.5 

eP 

25.7 

E 

0.5 

0.6 

iS 

26.2 

c 

0.3 

0.8 

iS 

27.5 

F. 

0.4 

0.6 

iS 

28.3 

N 

0.8 

0.4 

22 

1-15-56 

eP 

05-04-05.3 

N 

0.2 

pP 

05.8 

S 

0.3 

0.6 

eP 

05.8 

E 

0.2 

0.4 

iS 

07.0 

S 

0,7 

0.6 

iS 

07.4 

W 

0.7 

0.4 

23 

1-15-56 

eP 

05-03-05.9 

c 

0.2 

iP 

05.9 

S 

0.3 

0.4 

eP 

05.9 

E 

0.4 

- 

iS 

06.8 

N 

0.6 

0.4 

iS 

06.8 

W 

0.5 

0.4 

i 

07.2 

S 

0.9 

0.4 

6 


No. 

.uate 

Phase 

Time 

Direction 

Amplitude 

Period 

24 

11-23-55 

IP 

09-59-28.9 

S 

0-3 

0.4 

IP 

28.9 

W 

0.4 

0.4 

IS 

29.9 

S 

0.7 

0.4 

IS 

30.4 

W 

1.0 

0.4 

25 

10-15-55 

iP 

04-22-41.1 

d 

0.1 

0.4 

eP 

41.9 

W 

0.2 

0.4 

eP 

42.1 

S 

0.2 

- 

iS 

42.7 

E 

0.5 

0.2 

26 

10-15-55 

eP 

04-21-21.4 

E 

0.1 

0.2 

eP 

21.5 

c 

0.1 

0.4 

IP 

21.7 

N 

0.1 

0.3 

IS 

22.7 

E 

0.4 

0.2 

iS 

23.2 

N 

0.5 

0.6 

27 

9-17-55 

eP 

20-19-45.4 

d 

0.1 

0.2 

eP 

45.9 

N 

0.1 

- 

eP 

45.9 

W 

0.05 

- 

eS 

46.2 

d 

0.3 

0.2 

iS 

46.6 

S 

0.7 

0.6 

eS 

46.6 

E 

0.4 

0.6 

28 

9-2-55 

eP 

08-17-53-4 

c 

0.2 

0.6 

eP 

53-9 

E 

0.1 

0.4 

IP 

54.2 

S 

0-3 

0.6 

IP 

54.5 

c 

0.7 

0.5 

IP 

54.8 

N 

0.5 

l.o 

IP 

55.2 

E: 

0.7 

0.6 

IP 

55-2 

d 

0.8 

0.4 

iS 

56.1 

c 

0.7 

- 

iS 

56.1 

S 

l.l 

0.6 

iS 

56.5 

E 

l.o 

0.4 

29 

8-30-55 

eP 

02-19-26.7 

c 

0.1 

0.6 

IP 

26.7 

N 

0.3 

0.6 

eP 

26.7 

W 

0.05 

- 

iS 

29.5 

N 

0.8 

0.6 

eS 

30.2 

E 

0.6 

0.5 

eS 

30.2 

d 

0.3 

0.8 

i 

J2.4 

W 

0.9 

0.4 

30 

7-29-55 

IP 

12-48-15-3 

d 

0.8 

0.4 

IP 

15.9 

S 

0.7 

- 

IP 

15.9 

W 

1.2 

0.4 

7 


LiO  • 

Tate 

Phase 

Time 

Direction 

Amplitude 

Period 

51 

5-51-5: 

iP 

06-01-19.9 

d 

0-3 

- 

iP 

20.1 

S 

0.6 

0.4 

iP 

20.1 

E 

0.5 

0-5 

i3 

20.5 

W 

1.0 

0.6 

iS 

20.9 

N 

0.6 

0.4 

52 

2-11-55 

eP 

20-59-25*6 

d 

0.05 

- 

ep 

25.6 

S 

0.2 

0-3 

e? 

25.6 

E 

0.2 

- 

iS 

27.9 

W 

0.4 

- 

eS 

28.3 

N 

0.4 

0.3 

eS 

28.4 

c 

0.2 

0.6 

e 

29.1 

N 

0.5 

0.6 

55 

12-22-54 

e? 

00-35-28.6 

c 

0.05 

- 

eP 

29.1 

S 

0,1 

0.4 

iP 

29.4 

W 

0.1 

0,2 

iS 

29.8 

N 

0.4 

0.4 

iS 

30.5 

N 

o.4 

0.4 

eS 

30.6 

d 

0.2 

0.4 

IS 

31.0 

E 

0,1 

0.4 

i 

32.0 

S 

0.5 

0.6 

34 

9-14-54 

i? 

19-57-59.5 

d 

o.l 

- 

eP 

59*7 

N 

0.2 

- 

eP 

59-7 

E 

0.2 

0.6 

iS 

58-OO.9 

d 

0.2 

0.4 

iS 

01.4 

c 

0.3 

0.4 

eS 

01.6 

N 

0.3 

1.0 

eS 

C2.2 

E 

0.4 

1.0 

55 

6-30-54 

eP 

12-58-05,3 

E 

0.1 

- 

iP 

03.7 

d 

0.2 

- 

iP 

05.7 

a 

0,4 

- 

iS 

04.5 

c 

0.4 

- 

iS 

04. 5 

s 

1.2 

0.4 

iS 

o4.6 

w 

0.8 

0.4 

36 

8-27-53 

e? 

10-00-03,6 

s 

0.2 

0.6 

eP 

05.6 

w 

0.1 

0.6 

e? 

04. 3 

c 

0.05 

0.2 

eS 

05.7 

w 

0.2 

1.0 

iS 

06.3 

Q 

0.7 

0.5 

eS 

06.4 

c 

0.2 

0.6 

5T 

5-1-53 

iP 

35-37-48.7 

s 

0.2 

- 

iP 

48.8 

c 

0.2 

0.1 

iP 

48.9 

E 

0.2 

- 

iS 

49.4 

S 

0.9 

0.4 

iS 

49.5 

d 

0.8 

- 

iS 

49.6 

E 

0.8 

- 
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No. 


Date 


Phase 


Time 


Direction  Amplitude  Period 


38  3-13-53 


39  11-24-52 


4o  11-3-52 


41  8-7-52 


42  2-12-52 


43  4-7-51 


eP 

13-58-31.9 

s 

0.3 

1.0 

eP 

31-9 

w 

0.1 

0.2 

eP 

32.3 

c 

0.2 

0.2 

iS 

32-9 

W 

0.8 

0.6 

iS 

33-5 

c 

0.3 

0.3 

IS 

33.5 

0 

0.7 

0.4 

iP 

15-56-19.5 

d 

0.1 

0.2 

eP 

19.5 

E 

0.4 

1.0 

iP 

19.6 

N 

0.2 

0  4 

iS 

20.3 

N 

0.3 

0.4 

iS 

20.3 

W 

1.2 

0.8 

iS 

20.6 

d 

0.5 

0.2 

eP 

05-48-09.8 

N 

0.05 

0.8 

iP 

10.3 

E 

0.2 

0.4 

eP 

ll.l 

d 

0.1 

0.6 

IS 

11.0 

N 

0.2 

0.6 

iS 

11.4 

W 

0.4 

0.4 

PS 

11.9 

d 

0.8 

1.4 

eP 

18-00-02.9 

c 

0.05 

0.3 

iP 

03.1 

d 

0.1 

0.4 

iP 

04.0 

N 

0.2 

- 

iP 

04.5 

W 

0.5 

0.6 

iS 

04.1 

c 

0.4 

0.2 

iS 

04.7 

N 

0.3 

0.2 

iS 

05.4 

E 

0.6 

0.6 

eP 

15-34-57.2 

d 

0.05 

- 

eP 

57-9 

W 

0.3 

0.6 

eP 

58.3 

N 

0.2 

0.6 

eP 

58.7 

c 

0.3 

- 

iS 

35-00.4 

V 

0.7 

0.6 

iS 

01.4 

C* 

U 

1.2 

0.6 

i 

04.3 

W 

0.8 

1.4 

eP 

07-32- 59 '5 

N 

0.1 

0.8 

eP 

59-7 

W 

0.1 

0.6 

iP 

33-OO.O 

E 

0.2 

l.o 

eP 

00.2 

c 

0.05 

0.4 

iS 

oi.4 

N 

0.3 

0.8 

iS 

oi.4 

E 

0.5 

0.4 

eS 

m  r 

c 

0.1 

0.3 

i 

02.4 

W 

0.7 

0.6 
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No. 

Date 

Phase 

Time 

Direction 

Amplitude 

Period 

3-2-51 

eP 

17-20-19.6 

c 

0.1 

0.6 

iP 

20.7 

S 

0-3 

0.7 

eP 

20.8 

E 

0.3 

0.6 

eS 

21.9 

d 

- 

- 

iS 

22.0 

E 

0.5 

0.6 

iS 

22.4 

S 

0.5 

0.7 

i 

26.5 

E 

2.0 

0.5 

45 

1-29-49 

iP 

22-28-23.5 

S 

0.1 

0.3 

iP 

23.5 

W 

0.? 

0.2 

IP 

23.6 

c 

0.1 

0.2 

IS 

24.1 

N 

0.5 

0.4 

eS 

24  1 

E 

0.6 

- 

IS 

24.1 

c 

0.5 

0.2 

46 

12-8-48 

e 

20-37-35-5 

S 

0.05 

i 

46.8 

N 

0.2 

0.4 

i 

47.2 

W 

0.2 

0.4 

e 

47.4 

d 

0.1 

- 

47 

12-8-48 

eP 

16-18-20.9 

N 

0.5 

0.4 

eS 

22.9 

S 

0.4 

0.6 

eS 

22.9 

c 

0,05 

- 

-48 

10-7-48 

IP 

11-43-22.8 

W 

0,2 

_ 

IP 

22.9 

d 

0.2 

0.1 

IP 

23.0 

N 

0,2 

0.2 

49 

8-21-48 

eP 

15-26-41.2 

E 

0,1 

eP 

41.5 

d 

- 

- 

iP 

4l.6 

S 

0.5 

0,8 

iS 

43.0 

d 

- 

- 

eS 

43.2 

E 

0.5 

0.5 

iS 

3-6 

S 

0.5 

0.4 

i 

45.3 

W 

0.9 

0.3 

50 

5-23-48 

eP 

07-14-15-8 

c 

0.2 

0,5 

51 

4-20-48 

eP 

02-03-48.2 

d 

0.05 

• 

iP 

48.4 

E 

0.1 

0.5 

eS 

50.0 

E 

0.1 

0.4 

iS 

50.6 

S 

0.3 

0.4 

52 

4-20-48 

iP 

04-06-23.9 

W 

o.l 

0.4 

eP 

25.0 

S 

0.3 

0.5 

iP 

25.5 

w 

0.2 

0.7 

iP 

25.7 

d 

0.1 

0.2 

iS 

26.3 

S 

0.9 

0.6 

iS 

27.2 

E 

0.7 

0.6 

iS 

27.5 

c 

0.3 

0.4 

i 

29.7 

E 

1.0 

0.6 

10 


No. 

Date 

Phase 

Time 

Direction 

Amplitude 

Period 

53 

4-16-48 

eP 

17-05-06.0 

E 

0.5 

06.3 

c 

0.4 

0.4 

IP 

06.4 

N 

1.0 

0.2 

3.3  ANALYSIS  OF  FIRST  MOTION 

An  examination  of  the  above  data  reveal*  that  out  of  the  fifty-three 
events  there  were  forty-two  observation*  of  vertical  component  first  motion = 
!l%ese  were  distributed  as  eighteen  compression*  and  twenty-four  dilatations 
In  those  cases  where  no  vertical  component  was  recorded  but  where  horizontal 
components  indicated  an  azimuth  direction,  the  sense  of  the  motion  is  am¬ 
biguous  so  that  no  conclusion  could  be  drawn. 

Combining  the  direction  of  motion  and  the  compression-dilatation 
nature  of  the  first  motion,  there  were  thirty-four  cases  for  which  these 
characteristics  were  known.  Arranging  these  In  the  general  quadrants, 
northwest,  northeast,  southeast,  and  southwest,  as  indicated  by  the  metiers 
above  gives  a  picture  of  the  relative  frequency  of  dilatations  and  compress 

s ions  arriving  at  the  observatory  from  each  quadrant.  This  data  is  given. 

« 

in  the  following  table: 


PNTngWTBI . 

NW 

NE 

■-H 

sw 

Compression 

3 

4 

3 

3 

Dilatation 

5 

5 

_ 5 _ 

6 

The  table  shows  that  for  eight  earthquakes  lying  in  the  quadrant  north¬ 
west  of  the  observatory,  three  of  these  had  compressional  first  motion  and 
five  had  dilatational  first  motion.  Similar  results  are  seen  for  the  ether 
quadrants.  This  shows  that  for  this  type  of  low  Intensity  volcanic  quake 
in  the  nan-eruptive  period,  the  first  motion  is  not  uniformly  compressional. 


it 


In  fact,  it  shows  for  the  limited  number  of  quakes  which  occurred  a  slight 
statistical  preference  for  dilatational  first  motion. 

Distance  based  on  the  (S-P)  interval  and  a  longitudinal  wav°  velocity 
of  1*79  km/second  were  computed  and  plotted  on  a  graph  with  the  ot  'ervatory 
as  the  center.  No  apparent  relationship  between  distance,  direction,  and 
compression-dilatation  nature  of  the  first  motion  is  apparent.  The  results 
are  shown  in  Figure  1. 

J>.k  AMPLITUDE -DISTANCE  RELATIONSHIP 

With  the  distance  values  it  is  possible  to  inquire  whether  a  significant 
relationship  exists  between  the  observed  amplitude  of  motion  and  distance 
from  the  epicenter.  The  medial  test  was  used  to  determine  if  such  a  re¬ 
lationship  exists.  This  test  consists  of  plotting  the  amplitude  values  as 
ordinates  versus  the  distance  values  os  abscissae.  A  horizontal  medial  line 
is  drawn  so  that  half  the  points  fall  above  the  line  and  the  other  half 
below  the  line.  Next  a  vertical  medial  line  is  drawn  so  that  half  of  the 
points  lie  to  the  left  of  the  line  and  the  other  half  lie  to  the  right  of 
the  Line.  Theoretically  equal  numbers  of  points  should  fall  in  opposite 
quadrants  such  as  quadrants  one  and  three,  and  quadrants  two  and  four.  If 
no  significant  relationship  exists  between  the  parameters  the  points  will  be 
equally  disuributed  in  the  four  quadrants  or  approximately  so.  If  there  are 
too  few  points  in  any  quadrant  then  a  relationship  between  the  parameters  is 
inferred.  Points  falling  on  a  medial  line  are  disregarded. 

Amplitude -distance  values  for  28  events  were  taken  from  Table  I  using 
the  maximum  amplitude  value.  Referring  to  Quenoilles  table  of  significant 
levels  for  numbers  of  points  falling  in  any  quadrant,  the  number  depends  on 
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the  probability  level  which  is  assumed'  For  example,  using  28  points  and  a 
piobabxiit„  'it  significant  level  of  .01,  the  upper  and  lower  limits  are  11 
and  3  recpeci  ^ely .  Assuming  that  no  relationship  exists  between  the  pa¬ 


rameters,  the  probability  that  a  quadrant  will  have  as  many  as  11  or  as  few 
as  3  points,  the  upper  and  lower  limits,  is  less  than  .01.  A  portion  of 
Quenoilles  tables  are  reproduced  here  in  Table  2. 


TABLE  2  SIGNIFICANT  LEVELS  FOR  NUMBER  OF  POINTS  PALLING  IN  ANY  QUADRANT 


Lower 

Limit 

Upper  Limit 

Number  of  Points 

% 

If 

5$ 

if 

10  - 

n 

0 

0 

5 

5 

12  - 

13 

0 

0 

6 

6 

14  - 

15 

1 

0 

6 

7 

I 

1 6  - 

17 

1 

0 

7 

8 

18  - 

19 

1 

1 

8 

8 

20  - 

21 

2 

1 

8 

9 

22  - 

23 

2 

2 

9 

9 

2b  - 

25 

3 

2 

9 

10 

2  6  - 

27 

3 

2 

10 

11 

28  - 

29 

3 

3 

11 

11 

30  - 

31 

k 

3 

11 

12 

32  - 

33 

b 

3 

12 

13 

34  - 

35 

5 

b 

12 

13 

3b  - 

37 

5 

b 

13 

U' 

38  - 

39 

6 

5 

13 

±b 

bo  - 

bl 

6 

C 

14 

15 

42  - 

43 

6 

5 

15 

l6 

bb  - 

45 

7 

b 

15 

l6 

b6  - 

47 

7 

6 

lo 

17 

bQ  - 

49 

8 

Y 

1 

16 

IT 

50  - 

51 

8 

*"7 

I 

17 

18 

6o  - 

61 

10 

o 

y 

20 

21 

TO  - 

71 

12 

11 

23 

24 

8o  - 

8l 

15 

13 

25 

27 

90  - 

91 

17 

15 

28 

30 

100  - 

101 

19 

18 

31 

32 

no  - 

111 

21 

20 

34 

35 

14 


When  the  medial  test  was  applied  to  the  owentv-eight  points  the 
distribution,  reading  from  Quadrant  I  through  Quadrant  IV  was  5,8, 5 >6. 
This  indicates  no  significant  relationship  between  amplitudes  and  dis¬ 
tance.  The  medial  test  is  shown  in  Figure  2. 

Ratios  of  amplitude  over  period  were  computed  for  these  28  points 
and  are  given  here  in  Table  3- 


TABLE  3  AMPLITUDE  OVER  PERIOD  vs  DISTANCE,  VESUVIUS 
Event  No.  Amplitude/Period  Distance  (km) 


1 

O.63 

2.20 

9 

2.50 

1.96 

10 

2.66 

2  45 

13 

O.71 

10.58 

Ik 

3.00 

3-92 

15 

0.80 

7.60 

l6 

9.00 

6.86 

17 

1.00 

9.80 

20 

0.50 

1-37 

23 

2.25 

2.20 

2b 

1.75 

2.85 

26 

0.33 

3.19 

27 

1.16 

I.96 

29 

2.25 

6.86 

31 

2.00 

1.37 

32 

0.83 

5.68 

33 

0.83 

2.98 

38 

1.25 

3.83 

35 

2.00 

2.98 

36 

1.25 

5.15 

37 

2.25 

1.71 

38 

1.33 

2.85 

0.80 

1.96 

80 

0.57 

3-19 

kb 

4.00 

3.19 

85 

1.25 

1-37 

b9 

3.00 

8.81 

52 

1.60 

3.19 

'The  above  values  were  plotted  against  distance.  The  medial  test  was 
applied  to  determine  if  a  significant  relationship  might  now  exist  between 
"amplitude/period"  and  distance.  The  results  shown  in  Figure  3  give  a 
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v_n 


distribution  over  the  Quadrants  I  through  IV  of  6,8, 6,8.  While  this  is  a 
slight  improvement  it  cannot  be  considered  a  significant  relationship. 

The  above  test,  though  failing  to  show  significance,  does  show  a 
tendency  toward  significance  by  the  slightly  improved  quadrantal  distribu¬ 
tion  over  that  of  the  amplitude -distance  distribution.  In  the  hope  of 
improving  the  significance  and  as  a  further  check  on  the  amplitude -distance 
relationship,  the  trace  amplitudes  were  reduced  to  ground  motion  for  the  iG 
events.  This  was  done  by  dividing  the  trace  amplitude  by  the  magnification 
factor  e.t  the  period  of  the  observed  wave.  The  ground  motions  are  given  in 
Table  4. 


TABLE  4  GROUND  AMPLITUDE  vs  DISTANCE,  VESUVIUS 


Event  No. 


Ground  Amplitude  10 


-3 


mm 


1 

3.06 

9 

6.27 

10 

5.02 

13 

3.09 

14 

16.36 

15 

1.80 

l6 

33.53 

IT 

4.22 

20 

4.38 

23 

5.64 

2k 

4.39 

26 

2.07 

27 

4.34 

29 

6.34 

31 

5.02 

32 

3.11 

33 

3.11 

31 

4.03 

35 

5.67 

36 

3-14 

37 

5.64 

38 

6.24 

39 

7.65 

4o 

5.20 

kk 

14.3.4 

45 

3-14 

49 

6.42 

52 

5.66 
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The  ground  amplitudes  of  Table  k  were  plotted  against,  distance  and 
the  medial  test  again  was  applied  as  shown  in  Figure  k.  The  results  give 
a  7>7>7>7>  quadrantal  distribution  for  the  28  points  showing  the  complete 
lack  of  any  significant  relationship  between  amplitude  and  distance. 

3-5  MAGNITUDE  U IVESTIGATIDN 

A  somewhat  different  approach  was  attempted.  Variations  in  trace 
amplitude  are  observed  and  the  above  tests  indicate  that  Lhese  are  inde¬ 
pendent  of  distance.  A  number  of  factors  could  be  involved  to  explain 
qualitatively  at  least  why  these  variations  occur.  Factors  such  as  azimu¬ 
thal  variations,  differences  in  wave  path,  differences  in  instrumental 
response,  variations  with  change  in  source  location,  and  finally  variations 
in  energy  at  the  source.  Since  the  source  area  is  limited  and  the  distances 
are  all  small,  the  main  factor  would  seem  to  be  variation  in  energy  at  the 
source.  Assuming  all  other  factors  to  be  constant,  variations  in  energy 
at  the  source  will  become  manifest  as  variations  in  amplitude  of  the 
recorded  seismic  waves.  The  seismologist  deals  with  this  problem  by  com¬ 
puting  a  magnitude  for  the  earthquake  in  question. 

Richter,  in  dealing  with  local  California  quakes,  developed  an  empirical 
relationship  of  the  type: 

M  =  log  A  -  log  Aq  (l) 

Where  A  =  Amplitude  in  millimeters  with  which 
a  standard  torsion  seismometer  with 
characteristics  (T  =0.8,  V  =2800,  h  =  0.8) 
should  register  an  earthquake  of  magnitude 
zero. 

The  zero  magnitude  earthquake  is  further  defined  as  a  shock  that  would  pro¬ 
duce  a  trace  amplitude  of  me -thousandth  of  a  millimeter  on  a  standard  torsion 
seismometer  at  a  distance  of  100  kilometers . 
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Empirical  values  for  log  Aq  have  been  developed  by  Richter  for  Southern 
California.  These  values  are  regional  and  are  not  meant  to  be  applied  in 
other  arees.  However,  in  an  attempt  to  determine  approximately  the  magni¬ 
tude  involved  in  the  Vesuvian  quakes,  the  values  from  Richter's  table  were 
used. 

Figure  5  is  a  graph  of  log  amplitude  vs  log  distance  with  curves  for 
earthquake  magnitudes  1,  2,  and  3  computed  from  Richter's  values  for  log  Aq 
at  the  various  distances.  Since  the  magnification  of  the  Wiechert  seismo¬ 
graphs  at  Vesuvius  is  of  the  order  of  several  hundred,  the  trace  amplitude 
values  t'ere  multiplied  by  a  factor  of  ten  to  bring  the  magnification  up  to 
the  same  order  of  magnitude  as  that  of  the  standard  torsion  seismometer. 

The  response  curves  for  the  three  components  were  not  identical  and  this 
introduces  a  further  approximation.  The  results,  however,  are  considered 
to  be  of  the  right  crder  of  magnitude  and  indicate  that  the  magnitudes  are 
of  the  order  2.  This  represents  a  relatively  high  energy  level  for  this 
type  of  quake  and  indicates  that  seismic  activity  at  Mount  Vesuvius  is 
probably  of  a  much  higher  level  than  that  indicated  by  tue  limited  number 
of  quakes  recorded  from  19^+6  to  1961  by  the  low  magnification  Wiechert 
Seismographs.  High  sensitivity  seismographs  will  undoubtedly  reveal  many 
hundreds  of  lower  intensity  quakes.  A  magnitude  of  2  is  about  the  lower 
limit  of  perceptibility  and  probably  few  of  these  quakes  are  felt.  None 
of  the  quakes  were  recorded  anywhere  but  at  the  Vesuvian  Observatory. 

Several  examples  of  Vesuvian  records  are  reproduced  in  Figure  6.  The 
extremely  short  duration  and  high  frequency  character  of  such  shocks  is 
apparent  in  these  records. 
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FIGURE  G  VESUVIAN  SEISMOGRAM 
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4 

MESSINA 


4.1  THE  MESSINA  OBSERVATORY 

The  Seismological  Observatory  located  at  Messina,  Sicily  is  part  of  the 
Institute  of  Geophysics  and  Geodetics  of  the  University  of  Messina.  Professor 
Antonino  Girlanda  is  Director  of  the  Observatory  assisted  by  Dr.  Biagio 
Federico  and  a  staff  of  technicians. 

The  Messina  Observatory  is  approximately  90  kilometers  north  and  slight¬ 
ly  east  of  Mount  Etna.  It  is  also  approximately  300  kilometers  south  and 
slightly  east  of  Mount  Vesuvius.  It  is  unlikely  that  Vesuvian  quakes  would 
be  detected  at  this  distance  even  by  seismographs  of  intermediate  magnifica¬ 
tion  in  the  thousands.  Examination  of  Messina  records  failed  to  show  any 
activity  that  could  be  correlated  with  the  quakes  read  at  Vesuvius. 

The  greater  proximity  of  Mount  Etna  however,  gave  promise  that  seismic 
events  associated  with  it  might  be  recorded  at  Messina.  Mount  Etna  is  in  a 
state  of  almost  continuous  low  intensity  activity  centered  around  small  par¬ 
asitic  craters  near  the  summit.  No  significant  eruptions  had  occurred  in  the 
years  prior  to  the  examination  of  the  records . 

4.2  MOUNT  ETNA  -  ABORTIVE  ERUPTION 

In  December  of  1949,  however,  Mt.  Etna  appeared  to  be  preparing  for  a 
major  eruption.  A  number  of  quakes  were  recorded  at  Messina  which  appeared 
to  be  centered  near  the  volcano  and  an  eruption  appeared  imminent.  No 
eruption  occurred,  however,  and  this  has  been  referred  to  as  the  Abortive 
Eruption  of  December,  19^-9*  The  recordings  were  made  by  low  magnification 
Wiechert  instruments  and  indicate  fairly  high  source  energy  to  be  recorded 
at  the  distances  involved.  Readings  from  the  seismogram  are  given  in  Table  5, 
using  the  same  quantities  as  defined  for  Table  1. 
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TABLE  5  MESSINA  -  MOUNT  EEM  QUAKES 


No. 

Date 

Phase 

Time 

Direction 

Amplitule  Perioa 

^  (Sg-Pg) 

1 

12-2-49 

ePg 

04-31-38.5 

c 

0.1 

0.4 

66.1  km 

ePg 

58.5 

E 

0.1 

— 

eP 

32-00.1 

N 

0.1 

— 

log 

O6.3 

E 

0.3 

0.4 

iSg 

06.3 

S 

— 

— 

i 

18.1 

E 

0.7 

2.0 

2 

12-2-49 

iPg 

04-47-30.5 

c 

0.1 

1.2 

78.8  km 

ePg 

30.5 

S 

0.1 

1.4 

iP* 

31.7 

E 

0.2 

0.6 

eSg 

39-8 

E 

0.3 

0.6 

i 

48.8 

W 

0.5 

2.0 

3 

12-2-49 

ePg 

05-39-45.4 

c 

...... 

....  - 

88.9  km 

ePg 

46.2 

c 

0.1 

0.4 

iPg 

46.2 

E 

0.1 

0.4 

ePg 

46.9 

N 

0.1 

0.4 

ISg 

55-9 

W 

0.2 

0.6 

i 

40-08.7 

w 

0.8 

2.2 

4 

12-2-49 

iPg 

06-02-36.6 

c 

0.2 

0.6 

91.5  Ion 

ePg 

37.4 

E 

0.2 

0.6 

iPg 

37.7 

N 

0.2 

— 

i 

42.2 

W 

0.5 

— 

eSg 

47.4 

a 

0.2 

2.0 

iSg 

48.9 

N 

0.5 

1.4 

i 

54.7 

w 

2.0 

3.3 

5 

12-2-49 

ePg 

06-22-14.8 

c 

0.1 

-  -  - 

50.8  km 

iPn 

17.4 

a 

0.7 

1.6 

ePn 

17.8 

E 

0.3 

0.6 

iPn 

18.0 

s 

1.0 

1.2 

iSg 

20.8 

E 

1.1 

3.0 

i 

36.9 

a 

4.7 

3.2 

6 

12-2-49 

ePg 

06-43-41.3 

c 

0.1 

1.2 

59*3  km 

iSg 

48.3 

S 

C  .4 

1.5 

i 

44-03.7 

w 

0.6 

2.8 

7 

12-2-49 

e 

06-48-13.2 

w 

0.1 

1.4 

e 

23.9 

w 

0.4 

2.6 

8 

12-2-49 

ePg 

06-58-37.8 

E 

0.1 

1.0 

26.3  km 

eSg 

40.9 

s 

0.1 

1.0 

i 

44.5 

V 

0.3 

1.6 

9 

12-2-49 

e 

07-22-01.0 

c 

0.1 

1.0 

e 

01.4 

s 

0.1 

1.2 

i 

21.4 

E 

1.0 

2.4 
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No. 

Date 

Phase 

Time 

Direction  Amplitude 

Period 

£x  (Sg-Pg'1 

10 

.15 -2- 19 

e?g 

08-00-37. 1 

c 

0.2 

0.6 

59-3  km 

eP* 

37  .o 

N 

0.2 

1.2 

eP-:' 

38.2 

E 

0.1 

— 

iP* 

38.3 

c 

0.5 

1.0 

iSg 

77.1 

E 

0.8 

2.0 

iBg 

1+7.!+ 

E 

1.7 

2.7 

i 

56.8 

W 

5-0 

3-0 

11 

12 -2 -U9 

ePg 

08-15-26.6 

N 

0.1 

1.0 

iPg 

26.7 

d 

0.1 

— 

ePg 

26.9 

E 

0.2 

— 

i 

71 . 7 

c 

J 

0.7 

3-0 

1 

lb -17. 2 

“71 

ii 

7.0 

3-0 

The  Jeffrey's 

Pullen  Tables 

for  near 

earthquakes 

are  used 

at  Messina  and 

the  distances  listed  above  were  computed  from  these  trxles. 

Of  these  eleven  quakes  seven  began  with  initial  compressions!  phases, 
one  with  a  dilatation,  and  three  with  indeterminate  first  motion.  In  the 
case  of  the  dilatational  first  motion,  the  horizontal  component  motions 
indicate  a  source  to  the  Northeast  of  Messina  away  from  Mount  Etna.  For 
the  seven  compressional  first  motions,  three  have  horizontal  motions  indicat¬ 
ing  a  source  to  the  Southwest  in  the  direction  of  Mount  Etna.  The  distances 
are  88.9,  91- 5>  and  50.8  kilometers  respectively.  The  two  larger  distances 
would  pi  -i  the  source  near  the  cent''  r  of  Mount  Etna  while  the  shorter  dis  - 
tance  for  the  third  of  these  three  quakes  would  place  the  source  to  the 
Northeast  of  Mount  Etna.  Of  the  remaining  four  quakes  one  is  a  compression 
from  the  South  at  a  distance  66.1  kilometers  away,  again  falling  to  the  north 
of  Mount  Etna .  The  second  quake  was  from  the  north  of  Messina  and  the  two 
remaining  quakes  were  compressions  with  undetermined  direction. 

Those  data  are  very  sparse  and  any  conclusions  to  be  inferred  therefrom 
should  be  viewed  in  light  of  the  paucity  of  the  uata.  First  motion  from  the 
region  of  Mount  Etna  is  compressional  In  the  few  cases  observed.  Ot  er  quakes 
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which  should  probably  be  classified  as  regional  tectonic  quakes  originate 
in  this  general  area  off  the  flank  of  Mount  Etna.  These,  too,  displayed 
compressional  first  motion.  The  data  is  ;oo  scarce  however,  for  these 
observations  to  be  considered  as  more  than  observations. 

A  comparison  of  Messina  Observatory  records  for  quakes  originating 
near  the  center  of  Mount  Etna,  with  those  originating  to  the  north  of 
Mount  Etna  between  it  and  Messina,  show  a  rema  rki  ble  similarity  in 
character.  Seismograms  of  such  quakes  are  shown  in  Figure  7 • 

5 

MOUNT  ETNA 

5-1  EXPERIMENTAL  SEISMIC  INSTALLATION 

An  experimental  seismographic  installation  using  a  Willmore  three 
component  system  had  been  in  operation  on  Mount  Etna  just  prior  zc  this 
investigators  arrival  there.  Professor  Alfredo  Rittmann  directs  the 
Volcanological  Institute  of  Mount  Etna  at  the  University  of  Catania, 
Sicily.  The  Willmore  instruments  were  on  loan  to  the  Institute  through 
the  cooperation  of  Professor  J.  P.  Rothe  of  the  Institute  de  Physique 
du  Globe  in  Strasbourg,  France.  Through  the  courtesy  of  Professor 
Rittmann,  Professor  Rothe  and  Dr.  Haroun  Tazieff  of  the  Centre  National 
de  Volcanologie,  Belgium,  the  Willmore  records  were  made  available  for 
examination.  The  available  records,  however,  showed  very  limited  earth¬ 
quake  activity.  Further,  since  the  installation  was  experimental, 
instrumental  constants,  direction  f  motion  and  time  corrections  were 
usually  not  available. 
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FIGURE  7  MESSINA  SEISMOGRAMS 
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5.2  ANALYSIS  OF  EXPERIMENTAL  SEISMOGRAMS 

Six  quakes  were  read  with  S-P  values  ranging  from  1.2  to  1 . 5 
seconds  for  five  of  the  quakes  while  the  sixth  quake  had  an  S-P  of 
24.5  seconds.  This  latter  quake  can  be  dismissed  as  a  regional  tectonic 
quake  since  the  distance  is  or  the  o: ler  of  210  kilometers.  The  other 
five  quakes  ranged  in  distance  from  10.1  to  58*1  kilometers  and  represent 
hypocer.ters  located  outward  along  the  volcano's  flank.  No  conclusions 
regarding  first  motion  could  be  made  since  instrumental  response  was  not 
constant  due  to  the  experimental  nature  of  the  installation. 

The  most  prominent  characteristic  of  the  Mount  Etna  records  however, 

■was  the  persistent  tremor.  Tremor  is  usually  associated  with  an  eruptive 
period  and  Mount  Etna  is  in  a  state  of  almost  continuous  low  intensity 
activity.  The  tremor  had  a  remarkably  sinusoidal  character  showing  a 
prominent  beat  phenomena.  Readings  were  taken  of  the  direction  of  motion 
for  this  beat  phenomena  at  the  start  of  the  beat  at  five  stations  on  two 
different  days.  Direction  of  motion  was  constant  at  each  station  during 
the  period  when  readings  were  taken.  For  a  total  of  250  readings  there 
were  124  up  and  126  down  with  the  individual  station  readings  showing  similar 
distributions.  This  is  a  random  pattern  of  motion  and  indicates  that  the 
seismic  signal  is  probably  a  complex  mixture  of  motions  from  many  sources 
within  the  volcano,  and  not  a  simple  motion  resulting  from  a  pair  of 
frequencies.  A  V/illmore  Seismogram  from  Mount  Etna  is  shown  in  Figure  8. 
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6 


REGGIO  CALABRIA 

Records  examined  ai  Reggio  Calabria  through  the  courtesy  of 
Professor  Vittorio  3aror.e-Adesi,  Director  of  the  Observatory,  were  all 
registered  by  Wiechert  Seismographs.  The  seismograms  represented 
regional  and  local  seismic  activity  and  no  correlation  could  be 
established  with  volcanic  activity  at  Mount  Vesuvius  or  at  Mount  Etna . 
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7 


HAWAIIAN  VOLCANOES 


7-1  THE  HAWAIIAN  VOLCANO  OBSERVATORY 

The  Hawaiian  Volcano  Observatory  is  situaued  on  the  Island  of  Hawaii,  the 
largest  of  the  eight  major  islands  comprising  the  Hawaiian  group.  It  is  on 
the  northwest  rim  of  Kilauea  Caldera  approximately  500  feet  above  the  caldera 
floor.  Almost  due  south  of  the  Observatory  in  the  southwest  portion  of  the 
Caldera  is  Halemaumau,  the  fire  pit  or  central  vent  of  the  Kilauea  volcano. 
Almost  due  east  of  the  Observatory  is  a  smaller  crater,  Kilauea  Iki,  in  which 
an  eruption  occurred  in  November,  1959* 

At  that  time,  Dr.  Donald  H.  Richter  was  Scientist  in  Charge.  Through 
his  cooperation  and  assistance  and  with  permission  from  Dr.  Thomas  B.  Nolan, 
Director  of  The  United  States  Geological  Survey,  under  whose  auspices  the 
Hawaiian  Volcano  Observatory  is  operated,  seismic  data  from  the  Kilauea  Iki 
eruption  were  collected. 


7.2  THE  SEISMOGRAPH  NETWORK 

A  network  of  nine  stations  is  operated  on  the  Island  of  Hawaii  by  the 
Volcano  Observatory.  Directly  in  the  Kilauea  Caldera  is  the  North  Pit 
Station.  Just  outside  the  Caldera  within  a  distance  of  several  kilometers 
are  the  stations  Uwekahuna,  Outlet,  Desert,  and  Fauna  Loa.  With  the  ex¬ 
ceptions  of  Uwekahuna,  these  stations  record  by  telemetering  to  the  Observatory. 
At  distances  measured  in  tens  of  kilometers  are  the  stations  Hilo,  Naalehu, 

Pahoa  and  Kamuela.  These  stations  are  shown  on  the  map  in  Figure  9*  Station 
locations  and  instruments  are  given  in  Table  6. 
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FIGURE  9  SEISMIC  NETWORK,  HAWAIIAN  VOLCANO  OBSERVATORY 
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SABLE  6  H.  V.  0.  SEISMIC  NETWORK 


Station 

Coordinates 

Seismograph 

North  Pit 

19°  24.9' 

N 

HV0  -  2,  Z  Component 

155°  17-0' 

W 

Remote  recording 

h  =  1115 

m 

Uwe kahuna 

19c  25.4' 

N 

HVO  -  1,  Z  Component 

155°  17,8' 

W 

Springnether  Short  Period,  E,  Z  Components 

h  =  1240 

m 

Press -Ewing  Long  Period,  N,E,Z  Components 

Outlet 

19°  2J. 4' 

N 

HVO  -  2,  Z  Component 

155°  16. 9' 

W 

Remote  recording 

h  =  1080 

m 

Mauna  Loa 

19°  29.8' 

N 

HVO  -  2,  Z  Component 

155°  25.5' 

W 

Remote  recording 

h  =  2010 

m 

Desert 

19°  20.2* 

N 

HVO  -  2,  Z  Component 

155°  23.3' 

W 

Remote  recording 

h  =  815 

m 

Hilo 

19°  43.2' 

N 

HVO  -  1,  Z  Component 

155°  05-3' 

W 

Wood -Anders on,  N,  E  Components 

h  =  20 

m 

Pahoa 

19°  29,7' 

N 

Loucks-Omori,  N,  E  Components 

154°  56.8' 

W 

h  =  205 

m 

Naalehu 

19°  03.8' 

N 

Loucks-Omori,  N,  E  Components 

155°  35-2 

W 

h  =  205 

m 

Kamuela 

20°  01.3' 

N 

Loucks-Omori,  N,  E  Components 

155°  40.3' 

W 

h  =  815 

m 

Seismographs 

HVO-1,  electromagnetic,  T  =0.5  sec.,  maximum  magnification 
about  20,000 

HVO-2,  electromagnetic,  T  =0.8  sec.,  maximum  magnification 
about  20,000  0 
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Springnether,  electromagnetic,  T  =  1.3  sec.,  maximum  magnification 
about  2,000  0 

Wood -Anders on,  torsion  seismometer,  T  =  0.8  sec.,  maximum 
magnification  about  2,,800 

Loucks-Omori,  mechanical  seismograph,  T  =3*0  sec.,  maximum 
magnification  about  200 

7*3  Seismograph  Readings 

Kilauea  Iki  began  eruption  on  November  l4,  1959  at  about  8:00  P.M. 
local  time  or  20  h  GMCT.  Trior  to  the  eruption,  earthquake  activity  began 
and  increased  ’until  about  20  h  GMCT  when  the  eruption  occurred.  At  this 
time  earthquake  activity  as  such  ceased  and  was  replaced  by  volcanic  uremor, 
a  continuous  high  frequency  oscillation.  In  the  hours  preceding  the  erup¬ 
tion  many  hundreds  of  quakes  occurred,  so  many  that  frequently  they  became 
masked  and  unreadable  by  overriding  of  one  quake  onto  the  trace  of  another. 
Beginning  about  21  hr.  GMCT  on  November  13,  1959 >  readings  were  made  at  as 
many  stations  as  possible  up  to  the  time  of  the  eruption.  These  readings 
are  given  in  Table  7 


TABLE  7  SEISMOGRAPH  READINGS  -  KILAUEA  IKI  ERUPTION 


Station  Observations: 


U  =  Uwekahuna 
NP  =  North  Pit 
0  =  Outlet 
D  =  Desert 
ML  =  Mauna  Loa 
H  =  Hilo 


S  Phase  readings  are  tabulated  at  Uwekahuna  since  this  is  the  only  station 
with  high  sensitivity  horizontal  seismometers.  Readings  at  all  other 
stations  are  exclusively  P  readings  unless  otherwise  indicated. 


Initial  motion  is  indicated:  c  =  compression,  d  =  dilatation 
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Qua  ke _  t?«  ll«P.  0.  B.  M«L.  KT 


No." 

h  m 

T> 

s 

P 

P 

P 

P 

1 

21-39 

19.4  C 

20.4 

18.3  c 

19.6 

c 

21.9 

d 

2 

21-41 

05.2 

c 

3 

21-45 

10.6  d 

10.9 

c 

13.2 

c 

4 

21-50 

53-4 

d 

5 

22-00 

56.7 

c 

6 

22-05 

21.7 

d 

7 

22-09 

20.3  c 

20.9 

d 

21.7 

d 

22.1 

d 

8 

22-09 

25.9 

d 

9 

22-09 

42.2 

d 

10 

22-18 

27.9 

c 

li 

22-18 

46. 5  c 

47.1 

46.7 

d 

50.1 

c 

12 

22-42 

19.1 

d 

23.6 

c 

13 

23-23 

29.I 

c 

14 

23-35 

19.6 

c 

15 

23-39 

03.0 

- 

16 

23-50 

55.7  d 

56.5 

55.9 

d 

57-7 

c 

17 

oo-o4 

26.1  c 

18 

00-09 

56.6  c 

57-9 

49.0 

d 

50.6 

c 

19 

00-10 

57-5  C 

58.3 

20 

00-15 

16.3 

c 

21 

00-18 

38.7 

d 

22 

00-19 

26.2  c 

26.8 

28.2 

d 

28.6 

c 

23 

00-22 

12.9 

d 

24 

00-27 

31.7  d 

32.4 

25 

00-28 

06.4  d 

07.4 

08.3 

c 

08.8 

c 

26 

00-29 

52.3  d 

53-2 

27 

00-30 

lb.4 

d 

28 

00-32 

32.5  d 

33.5 

32.1 

c 

35.3 

d 

29 

00-35 

45.0  d 

45.9 

44.9 

c 

48.4 

d 

30 

OO-36 

08.8  d 

10.2 

09.4 

c 

11.4 

c 

31 

00-37 

51.2  d 

51.6 

56.5 

c 

32 

00-39 

45.2  c 

46.3 

48.9 

d 

33 

00-39 

58.5  d 

59-3 

34 

00-42 

06.7  d 

07.8 

06.6 

c 

09.7 

d 

35 

00-43 

25-4  c 

26.0 

36 

00-44 

23.6 

c 

37 

00-45 

41.6 

d 

38 

00-46 

09.6  c 

10.7 

11.5 

c 

39 

00-47 

15.5  d 

16.2 

4o 

00-47 

50.9  c 

51.4 

41 

00-48 

50.9  d 

51.4 

53.4 

d 

53-6 

c 

42 

00-49 

53.3  d 

54.3 

43 

00-51 

33.1  c 

33-6 

35.2 

d 

35.6 

d 

44 

00-52 

44.0  c 

45.4 

47.0 

c 

45 

00-54 

02.0  d 

02.9 

02.1 

c 

05.9 

c 

46 

00-54 

43.8 

c 

47 

00-54 

55.4  c 

56.3 

54.8 

c  02.4  c 

48 

00-55 

56.0  c 

57.1 

56.3 

d 

58.5 

c 
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Quake 

u. 

N,P. 

0 

D 

M.L. 

H. 

6 

h  m 

p 

s 

P 

P 

P 

P 

P 

49 

OO-58 

05.0  d 

05-7 

50 

00-59 

09.5  d 

51 

01-00 

49.2  d 

5^.0 

52 

01-01 

49.2  d 

50.0 

51.5 

c 

55 

01-03 

13.4  d 

14.4 

12.4 

c 

54 

01-03 

50.0 

c 

55 

01-05 

44.4 

- 

56 

01-11 

23-5  d 

24.6 

57 

01-11 

59.9 

d 

58 

01-12 

25-2  d 

26.1 

26.9 

d 

28.5 

d 

59 

01-14 

38.5  d 

39*2 

40.5 

d 

41.2  c 

60 

01-21 

58.4  d 

59*3 

58.9  c 

58.3 

c 

61 

01-23 

29.6 

d 

30.3  d 

62 

01-28 

12.6  d 

13.4 

11.9  0 

65 

01-28 

39-7  d 

40.5 

39-2  C 

64 

01-31 

08.0  C 

08.8 

07.4  c 

65 

01-31 

44.0  d 

45.0 

45.1  c 

66 

01-33 

15.4  d 

6? 

01-41 

04.1  c 

0^.1 

06.9 

d 

07.0  c 

68 

01-41 

30.0  c 

31.3 

29.4  c 

69 

01-42 

3-3.0  d. 

38.9 

37.4  c 

38.5 

d 

70 

01-43 

52.8  c 

33-8 

32.3  C 

32.8 

c 

35.7  c 

71 

01-44 

17-6  c 

17.9 

c 

21.6  c 

72 

01-45 

24.3  d 

25-0 

23.5  d 

23.7 

c 

73 

01-49 

10.3  d 

11.2  c 

74 

01-49 

56.9  c 

57.6 

55  9  c 

56.5 

c 

58.7 

c 

59-7  c 

75 

0.1-53 

47.4 

c 

76 

01-, 

56.4  c 

77 

01-55 

36.7  d 

37.4 

35-9  c 

36.4 

r* 

78 

01-56 

38.3  d 

39-1 

79 

02-00 

53.4 

c 

80 

02-01 

26.2 

c 

81 

02-05 

22.8  c 

23.6 

22.4  d 

23.1 

c 

24.9 

d 

25.2  d 

29  •  5  c 

82 

0.1 -06 

47.5 

d 

83 

02-09 

16.1 

c 

84 

02-12 

59-6  d 

00.3 

59.0  c 

59^ 

c 

01.6 

d 

02.4  d 

85 

02-16 

42.7 

d 

86 

02-18 

54.9  c 

56.I 

87 

02-22 

41.4 

c 

88 

02-24 

35.3 

c 

89 

02-26 

43.4  c 

44.3 

43.1  d 

43.7 

d 

46.4  c 

90 

02-28 

38.8  d 

39.7 

38.3  c 

91 

02-33 

47.1  c 

48.6 

46.4  d 

49.2 

c 

49.2  c 

vn 

CO 

P* 

92 

02-34 

31.4  c 

93 

02-36 

22.2  c 

94 

02-36 

49.2  c 

95 

02-37 

59-1  c 

00.1 

56.8  c 

59-4 

c 

00.6 

d 

01.6  c 

96 

02-39 

19-3  c 

20.4 

18.0  c 

37 
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No. 

h  m 

F 

P 

p 

P 

p 

P 

97 

02-41 

25.2  d 

26.1 

24.5  c 

25-3 

c 

27.4  d 

28.1  c 

98 

02-1:2 

03.4  c 

90 

02-44 

07.6 

c 

100 

02-44 

44.5  d 

101 

02-48 

17.0  c 

17.4 

d 

19-3  1 

19.5  c 

102 

02-49 

11.4 

c 

103 

02-37 

16.7  c 

16.7 

c 

18.9  d 

101+ 

03-00 

22.7  c 

23.8 

22.6  d 

23.0 

c 

26.7  c 

25.6  d 

105 

03-01 

12.7  a 

13.5 

12-3  c 

106 

03-03 

26.2  d 

107 

03-10 

19.9 

d 

108 

03-12 

00.2 

c 

109 

03-12 

52.0  d 

110 

03-13 

15.6  d 

16.5 

15.2  c 

15.7 

c 

111 

03-15 

36.2  d 

37-5 

35-7  1 

36.1 

c 

38.0  c 

38.2  d 

112 

03-17 

36.5  c 

37-3 

35-9  1 

38.9  c 

39.0  d 

43.0  d 

113 

03-21 

39-2  c 

n4 

03-22 

09.6  d 

113 

03-22 

24.6  d 

116 

03  23 

24.6  d 

25-1 

23.9  c 

24.6 

c 

26.4  d 

28.1  c 

117 

03-24 

I4.7  c 

16 . 0 

l6.0  c 

17.1  d 

21.6  d 

118 

03-25 

51.5  c 

53 .8 

c 

119 

03-27 

45.9 

c 

120 

03-28 

30-8  c 

31.4 

30.2  c 

33.2  c 

33.3  d 

370  c 

121 

05-31 

O6.6  d 

08.0 

06.4  d 

08.5  c 

12.6  c 

x22 

03-31 

37-9  c 

123 

03-35 

27.0  c 

121* 

0;  38 

50.8  d 

53-3 

50.4  c 

53.2  d 

5^.6  c 

59.0  d 

123 

03-40 

53-1  c 

54.1 

52.5  c 

54.4  c 

55.2  a 

58.8  d 

126 

03-43 

50.4  c 

51.2 

50.1  d 

127 

03-48 

29.9  d 

30.3  0 

35.7  c 

728 

03-51 

55-8 

a 

57.1  c 

129 

03-54 

23-9 

c 

130 

03-54 

53.6 

a 

131 

03-55 

52.5 

d 

132 

03-56 

54.2 

C 

13.3 

03-57 

13.1 

c 

1^4 

04-01 

3.I.9 

d 

135 

04-02 

49.0  c 

50.0 

48.5  c 

50.2  d 

51.6  d 

136 

04-04 

J16.O  c 

3  37 

04-05 

05.0  d 

138 

04-05 

24.6  c 

139 

04-07 

O8.3  c 

140 

04-08 

22.1 

c 

142 

04-1; 

00.6  c 

02.0  c 

142 

04-10 

22.6 

c 

38 


Quake 

U 

N.P. 

0. 

D 

M.L. 

H. 

No. 

h  m 

p 

S 

P 

P 

P 

p 

P 

143 

04-11 

42.5 

d 

144 

04-15 

25.1 

d 

26.0 

25-6 

d  27.2 

d 

27.7  C 

145 

04-17 

30.2 

c 

146 

04-13 

46.0 

c 

47.3 

45.7 

d 

47.0 

d 

48.8  c 

l4ff 

04-23 

20.8 

c 

22.1 

20.7 

d 

23.4 

d 

23.5  a 

148 

04-24 

52.5 

d 

53-6 

52.2 

d 

52.2 

c 

54.8  d 

'49 

04-25 

32.9 

c 

150 

04-2o 

45.6 

c 

46.6 

45.2 

d 

45.4 

c 

47.1 

c 

48.1  d 

51.5 

c 

151 

o4-28 

16.7 

c 

152 

04-32 

18.4 

d 

153 

o4-34 

18.4 

i 

154 

04-34 

39.7 

c 

155 

04-3J,' 

54.3 

c 

55.8 

54.3 

d 

54.4 

d 

156 

04-35 

59.9 

c 

157 

04-38 

14.9 

c 

158 

04-39 

16.5 

c 

17-7 

16.9 

c 

18.9 

d 

19.1  d 

23-2 

d 

159 

04-43 

02.1 

c 

03.6 

02.4 

c 

04.3 

d 

04.7  d 

160 

04-44 

01.4 

d 

01.5 

c 

161 

04-44 

26.0 

d 

26.8 

25*7 

c 

25.9 

c 

28.9 

d 

28.1  d 

162 

04-45 

18.4 

d 

163 

04-48 

38.9 

c 

164 

04-48 

55-4 

d 

165 

04-49 

23.4 

d 

25*3 

c 

166 

04-51 

17.2 

d 

21.3 

c 

167 

04-52 

25.9 

d 

28.5 

c 

168 

04-52 

32.6 

c 

33-8 

32.5 

d 

35-1 

c 

35.1  d 

38.9 

c 

169 

04-56 

07.6 

c 

08.1 

d 

170 

04-59 

56.5 

c 

171 

05-00 

33-7 

d 

172 

05-01 

08.9 

c 

173 

05-04 

39-2 

d 

4o.8 

39-4 

d 

174 

05-05 

03.4 

c 

o4.6 

04.2 

d 

05.1 

c 

05.5  d 

09.2 

c 

175 

05-06 

38.5 

c 

176 

05-07 

22.4 

c 

177 

05-07 

58.9 

d 

59-7 

58.6 

c 

173 

05-10 

06.5 

c 

07.6 

179 

05-14 

16.3 

d 

17.6 

19.3  c 

180 

05-15 

40.3 

c 

41.2 

181 

05-18 

36.1 

c 

36.9 

37-9 

c 

38.2  d 

182 

05-22 

370 

d 

38.6 

d 

39-7 

c 

39-9  a 

183 

05-24 

08.4 

c 

184 

05-24 

33.9 

d 

185 

05-26 

05.8 

d 

13.4 

c 

09.3  d 

17.? 

d 

186 

05-28 

25.0 

c 

187 

05-50 

15.8 

c 

19.7  d 

188 

05  .0 

27.1 

c 

189 

05-31 

30.4 

c 
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u 

N.P. 
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D. 
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H. 

No. 

n  m 

P 

S 

P 

P 

P 

P 

P 

190 

05-33 

17.6 

c 

18.7 

20.5  d 

191 

05-33 

43.6  c 

102 

05-3^ 

09.9  c 

193 

03-38 

36.9  c 

194 

05-42 

06.6  d 

195 

05-43 

41.7 

c 

45.6 

196 

05-4? 

05.7 

d 

06.7 

07-2  d 

07-8  c 

11.2  c 

197 

05-48 

23.1  c 

198 

05-49 

17-2  d 

199 

05-50 

08.2  c 

200 

05-56 

46.0  c 

201 

05-58 

12.9  c 

202 

06-01 

11.7  d 

203 

06-03 

05.3  c 

204 

06-03 

56.2 

d 

57.1 

59-1  d 

205 

Ob-06 

39-7  c 

206 

06-09 

51.8 

c 

54.2  d 

54.0  d 

58.3  d 

207 

06-12 

28.4 

c 

29-2 

30.9  d 

35.2  c 

208 

06-14 

08.0  d 

209 

06-19 

55-7  d 

210 

06-39 

29.4  d 

211 

07-35 

17.4 

c 

18.9 

20.2  c 

212 

0S-20 

51.8 

c. 

53-4 

54.9  c 

213 

08-33 

13-7 

c 

15.2 

16.6  d 

214 

08-41 

03.3 

d 

04.3 

07.2  c 

215 

10-16 

53-2 

c 

54.5 

56.O  c 

216 

17-59 

47.7 

d 

48.4 

The  acti . ity  at  the  North  Pit  Station  far  exceeded  that  of  any  other 
station.  From  22h  GMCT  until  Clh  2um  GMCT  nothing  could  he  read  because  of 
the  number  and  size  of  the  quakes.  After  Olh  20  m  GMCT  the  magnification  of 
the  instrument  was  reduced  by  a  factor  of  ten  from  20,000  to  2,000  after 
which  the  records  became  readable. 

A  simple  numerical  count  of  first  motion  dilatations  and  compressions 
at  each  station  is  given  in  Table  8.  The  distributions  for  Uwekahuna,  Maun a 
Loa,  and  Hilo  are  random.  At  North  Pit  and  Outlet  there  is  a  significant 
preponderance  nf  compress ional  motion.  At  Desert  a  veak  preference  for  dila- 
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TABLE  8  FIRST  MOTION  COMPRESSION  VS  DILATATION 


Station 

U 

NP 

0 

D 

ML 

H 

Compression 

31 

55 

01 

26 

39 

9 

Dilatation 

50 

„  - 

5^ 

29 

38 

8 

tational  motion  is  apparent.  This  possibly  indicates  a  critical  zone  in 
which  motion  may  originate  with  a  preferential  direction. 

A  group  of  quakes,  sixty-four  in  number,  was  selected  on  fche  basis 
of  being  recorded  at  three  or  more  stations.  First  motion  compressions 
and  dilatations  were  tabulated  for  this  group  giving  the  results  shown 
in  Table  9* 


TABLE  9  SELECTED  QUAKES  -  FIRST  MOTION  DISTRIBUTION 


No .  of 

No.  of 

Station 

No.  of  Quakes 

Compressions 

Dilatations 

Uwekahuna 

58 

30 

28 

North  Pit 

36 

19 

17 

Outlet 

39 

27 

12 

Desert 

hi 

16 

25 

Mauna  Loa 

5^ 

26 

28 

Hilo 

18 

10 

8 

This  distribution  shows  a  preferential  pattern  at  Outlet  and  Desert 
and  a  random  pattern  at  Uwekahuna,  North  Pit,  Mauna  Loa,  and  Hilo.  The 
preferential  pattern  at  Outlet  shows  a  dominance  of  compress ional  first 
motion,  while  at  Desert  there  is  a  dominance  of  dilatational  first  motion. 

Regrouping  these  data  into  two  sets,  one  with  first  motion  exclusive¬ 
ly  compressional  at  Uwekahuna  and  the  other  with  first  motion  exclusively 
dilatational  at  Uwekahuna,  the  first  motion  distributions  at  the  remaining 
stations  were  reexamined.  These  groupings  are  shown  in  Tables  10  and  11. 
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TABLE  10 


SELECTED  QUAKES  -  UWEKAHUNA  FIRST  MOTION  COMPRESSIONAL 


Station 

No.  of  Quakes 

No.  of 
Compressions 

No.  of 
Dilatations 

Uwekahuna 

30 

30 

0 

North  Pit 

l6 

5 

11 

Outlet 

1 6 

10 

6 

Desert 

21 

9 

12 

Mauna  Loa 

29 

12 

17 

Hilo 

11 

6 

5 

TABLE  11 

SELECTED  QUAKES  - 

UWEKAHUNA  FIRST  MOTION  DILATATIONAL 

No.  of 

No.  of 

Station 

No.  of  Quakes 

Compressions 

Dilatations 

Uwekahuna 

28 

0 

28 

North  Pit 

15 

10 

5 

Outlet 

20 

lo 

4 

Desert 

15 

5 

10 

Mauna  Loa 

21 

12 

9 

Hilo 

k 

2 

2 

Table  10  for  compress ional  motion  at  Uwekahuna  shows  a  significant 
dominance  of  dilatational  motion  at  North  Pit  and  of  compressional  motion 
at  Outlet.  Mauna  Loa  shows  a  weaker  dominance  of  dilatational  motion. 

Table  11  for  dilatational  motion  at  Uwekahuna  now  shows  a  reversal  at 
North  Pit  with  a  dominance  of  compressional  motion.  At  Outlet,  however, 
the  compressional  motion  dominance  is  very  significantly  strengthened. 

Desert  which  showed  only  a  slight  dominance  of  dilatational  motion  in 
Table  10  is  also  strengthened. 

These  observations  lend  strength  to  the  concept  of  a  critical  zone 
in  which  seismic  motion  may  originate  with  preferential  first  motion  patterns. 

The  relative  positions  of  these  quakes  can  be  somewhat  indicated  by  the 
S-P  values.  These  were  computed  for  Uwekahuna  for  the  57  quakes  in  the 
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selected  group.  The  S-P  values  are  given  in  Table  12. 


TABLE  12  S-P  VALUES  AT  UWEKAHUNA 


S-P 


No.  of  Quakes 


0.1 

0.2 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1-5 

1.4 
1-5 

2.5 


1 

1 

3 

4 

5 
5 

11 

9 

4 

3 

5 
2 

3 

1 


More  than  half  of  the  quakes  fall  in  the  S-P  range  0.8  to  1.1 
which  would  produce  a  distance  variation  of  approximately  one  kil¬ 
ometer.  Hence,  the  quakes  should  to.  i  to  cluster  together.  Varia¬ 
tions  in  S-P,  however,  may  indicate  increase  In  depth  rather  than 
an  increase  in  horizontal  distance.  In  order  to  evaluate  the  real 
significance  of  first  motion,  distance  and  position  of  the  station 
relative  to  the  earthquake  location  must  be  known. 


7,4  THE  EPICENTER  PROGRAM 

In  locating  epicenters  of  local  earthquakes  one  must  know 
something  about  the  crustal  structure  and  seismic  velocities  of 
the  area.  Scientists  at  the  Volcano  Observatory  have  developed 
a  crustal  model  with  appropriate  P  and  S  wave  velocities  that  fit 
the  seismic  observations  of  the  Hawaiian  area.  This  structural 
picture  is  sho^m  in  Figure  10.  Using  this  structure  and  assuming 
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FIGURE  10  CRUSTAL  STRUCTURE,  ILAlW.II 


zero  focal  depth,  travel  time  equations  have  been  computed  as  follow 


TP1  =  0.26  A 

TS  =  o.4iA 

2 

TP  =  1.0  +  0.20  A 

2 

TS  =  1.6  +  0.32  A 

2 

3 

TP  =  4.4  +  0.12  ^ 

3 

TS  =  6.8  +  0.21A 

3 

4 

Travel  time  curves  based  on  these  equations  are  shown  in  Figure  11. 
The  P  phase  does  not  appear  before  A  =  7*7  km  and  P  not  before 

2  3 

17.65  km.  Similarly,  S  does  not  appear  before^  =  7-6  km  and  S 

2  4 

not  before  20.2  km.  Travel  time  tables  for  zero  focal  depth  are 

given  in  Table  13. 


TABLE  13  TRAVEL  TIMES  FOR  ZERO  FOCAL  DEPTH 
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0.4 

.102 

.166 

0.5 

.128 

.207 
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1.04 
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FIGURE  11  TRAVEL  TIME  CURVES,  ZERO  FOCAL 
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Travel  times  were  also  computed  for  focal  depths  of  3>  4,  8,  12. 5, 
and  52.5  kilometers  and  these  values  are  given  in  Table  14.  Travel 
time  curves  based  on  these  computed  values  are  shown  in  Figure  12,  for 
the  P-phase  arrivals. 
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TABLE  l4  TRAVEL  TIMES  FOR  FOCAL  DEPTH 
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9-27 

13.7 

3.51 

5.71 

3.60 

5.83 

3.24 

5.25 

3.46 

5.59 

3.94 

7.02 

4.49 

7.42 

5.54 

9.28 

33.8 

3.54 

5-75 

3.62 

5.89 

3.26 

5.29 

3.48 

5.62 

3.96 

7.05 

4.50 

7.44 

5.55 

9.29 

13.9 

306 

5.79 

3.64 

5.91 

3.28 

5.32 

3.50 

5.65 

3.98 

7.08 

4.50 

7.45 

5.55 

9.30 

14.0 

3-59 

5.83 

3.67 

5.95 

3.30 

5-35 

3.52 

5.68 

3-99 

7.10 

4.51 

7.47 

5.56 

9.30 

14.1 

3.62 

5.88 

3-70 

5.99 

3-32 

5.38 

3.54 

5.71 

4.00 

7.13 

4.52 

7.49 

5.57 

9.32 

14.2 

3.64 

5.92 

3.72 

6.03 

3.34 

5.4i 

3.55 

5.74 

4.02 

7.16 

4.52 

7.50 

5.57 

9.32 
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5 


k 


S 


8 


h  0 

A  in  Ion  _ P  _S 

14. 3  3.67  5-96 

14.4  3.69  6.00 

14.5  3*72  6.04 

14.6  3.74  6.08 

14.7  3-77  6.12 

14.8  3-79  6.17 

14.9  3.82  6.21 
15.O  3.85  6.25 

15.1  3.87  6.29 

15.2  3.90  6.33 

15.3  3-92  6.38 
l$.4  3.95  6.42 

15.5  3.97  6.46 

15.6  4.00  6.50 

15.7  4.03  6.34 

15.8  4.05  6.58 

15.9  4.o8  6.62 
16.0  4.10  6.67 

16.1  4.13  6.71 

16.2  4.15  6.75 

16.3  4.18  6.79 

16.4  4.20  6.83 

16.5  4.23  6.88 

16.6  4.26  6.92 

16.7  4.28  6.96 

16.8  4.31  7-00 

16.9  4.33  7.04 
17.0  4.36  7.08 

17.1  4.38  7-12 

17.2  4.4l  7-17 

17.3  4.44  7.21 

17.4  4.46  7.25 

17.5  4.49  7.29 

17.6  4.51  7.33 

17.7  4.54  7.38 

17.8  4.56  7-42 

17.9  4.59  7.46 
18.0  4.62  7-30 

18.1  4.64  7.54 

18.2  4.67  7.58 

18.3  4.69  7.62 

18.4  4.72  7.67 

18.5  4.74  7.71 
3.8.6  4.77  7.75 

18.7  4.79  7-79 

18.8  4.82  7.83 

18.9  4.85  7-88 
19.O  4.87  7-92 


_P  _S  JP 

3.74  6.08  3.36  5.45 

3.77  6.12  3-38  5.48 

3.80  6.16  3.40  5.51 

3.82  6.20  3-43  5-54 
3.84  6.24  3.45  5.57 

3.87  6.29  3.47  5.61 

3.90  6.32  3.49  5*64 

3- 92  6.37  3.51  5.67 
3.94  6.41  3.53  5.70 
3.97  6.45  3-55  5.73 
4.00  6.49  3*57  >*77 
4.02  6.53  3.59  5  80 
4.o4  6.58  3.60  5.83 
4.07  6.62  3.62  5.86 
4.10  6.66  3.64  5.89 
4.12  6.70  3.66  5.93 
4.14  6.74  3.68  5.96 
4.17  6.78  3.70  5-99 
4.20  6.82  3.72  6.02 
4.22  6.86  3.74  6.05 
4.25  6.90  3.76  6.09 
4.27  6.94  3.78  6.12 
4.30  6.98  3.80  6.15 
4.33  7.03  3.82  6.18 

4- 35  7.07  3-84  6.21 
4.38  7-11  3-86  6.24 
4.4o  7-15  3-88  6.27 
4.43  7-19  3.90  6.30 
4.46  7.23  3.92  6.34 
4.48  7.27  3.94  6.37 
4.50  7-31  3.96  6.4o 

4.53  7.35  3.98  6.43 
4.56  7-40  4.00  6.46 
4.58  7-44  4.02  6.49 
4.6o  7.48  4.04  6.53 
4.63  7.52  4.06  6.56 
4.66  7.56  4.08  6.56 
4.68  7.60  4.10  6.56 

4.70  7.64  4.12  6.57 
4.73  7.68  4.14  6.57 
4.76  7.72  4.16  6.57 

4.78  7.76  4.18  6.57 

4.80  7.80  4.20  6.57 

4.83  7.85  4.22  6.58 
4.86  7.89  4.24  6.58 

4.88  7-93  4.26  6.58 

4.90  7.97  4.28  6.58 
4.93  8.01  4.30  6.59 


12.5 

PSPS 


3.57 

5-77 

4.04 

7.18 

3.59 

5-80 

4.05 

7-21 

3.61 

5.83 

4.06 

7.24 

3.63 

5.87 

4.08 

7.27 

3.65 

5.90 

4.10 

7.30 

3.67 

5.93 

4.11 

7.32 

3.69 

5*95 

4.12 

7.35 

3.70 

5.99 

4.14 

7.38 

3.72 

6.02 

4.16 

7.41 

3.74 

6.05 

4.17 

7.43 

3-76 

6.08 

4.19 

7.46 

3*78 

6.11 

4.20 

7.49 

3.79 

6.14 

4.22 

7.51 

3.81 

6.18 

4.24 

7.54 

3.83 

6.21 

4.25 

7.57 

3.85 

6.24 

4.26 

7.60 

3.87 

6.27 

4.28 

7.63 

3.88 

6.30 

4.30 

7.65 

3.90 

6.33 

4.32 

7.68 

3.92 

6.36 

4.34 

7.71 

3.94 

6.39 

4.35 

7.74 

3.96 

6.42 

4.37 

7-77 

3.98 

6.46 

4-39 

7.80 

4.00 

6.49 

4.41 

7.84 

4.01 

6.52 

4.43 

7.87 

4.03 

6.55 

4.44 

7-90 

4.05 

6.58 

4.46 

7-93 

4.07 

6.61 

4.48 

7.96 

4.09 

6.64 

4.50 

7.99 

4.11 

6.67 

4.51 

8.02 

4.13 

6.70 

4.53 

8.06 

4.15 

6.75 

4.54 

8.09 

4.16 

6.77 

4.56 

8.12 

4.18 

6.80 

4.58 

8.15 

4.20 

6.83 

4.59 

8.18 

4.22 

6.86 

4.6l 

8.22 

4.24 

6.89 

4.62 

8.25 

4.26 

6.92 

4.64 

8.28 

4.28 

6.95 

4.66 

8.31 

4.30 

6.98 

4.67 

8.34 

4.32 

7.01 

4.69 

8.37 

4.34 

7.04 

4.71 

8.40 

4.36 

7.08 

4.72 

8.43 

4.38 

7.11 

4.74 

8.47 

4540 

7.14 

4.76 

8.50 

4.42 

7.17 

4.78 

8.53 

4.44 

7.20 

4.79 

8.56 

4.46 

7.23 

4.81 

8.59 

22.5  32.5 

_p  JL  _p  _s 

4.53  7.32  5.58  9.32 

4.54  7.53  5.58  9.32 

4.55  7.55  5-59  9.32 

4.55  7. 57  5.60  9.32 

4.56  7.58  5.60  9-33 

4.57  7.60  5.61  9.35 

4.57  7.6l  5.61  9.33 

4.58  7.63  5.62  S-33 

4.59  7.65  5.63  9-33 

4.60  7.66  5.63  9.33 

4.60  7.68  5-64  9.34 

4.61  7.69  5.64  9.34 

4.62  7.71  5.65  9-34 

4.63  7-73  5.66  904 

4.64  7.74  5.66  9.34 

4.64  7.76  5.67  9-35 

4.65  7.77  5.67  9.35 

4.66  7.79  5.68  9.35 

4.67  7-79  5.69  9-35 

4.68  7.80  5.69  9.35 

4.68  7-81  5.70  9-35 

4.69  7.81  5.70  9-35 
4.7c  7-81  5.71  9.36 

4.71  7.82  5.72  9-36 

4.72  7-86  5-72  9.36 

4.72  7.89  5.73  9.36 

4.73  7-92  5.73  9-36 

4.74  7.96  5.74  9.56 

4.75  7.98  5-75  9.36 

4.76  7-99  5.75  9.36 

4.76  8.01  5-75  9.37 

4.77  8.03  5.76  9.37 

4.78  8.04  5.77  9.37 

4.79  8.06  5-77  907 

4.80  8.08  5.77  9.37 

4.80  8.10  5.78  9.38 

4.81  8.11  5-79  9.38 

4.82  8.13  5.79  9.38 

4.83  8.15  5.80  9.38 

4.84  8.16  5.80  9.38 

4.84  8.18  5.81  9.38 

4.85  8.20  5.81  9.38 

4.86  8.21  5.82  9-38 

4.87  8.23  5.83  9.39 

4.88  8.25  5.83  909 
“.88  8.27  5-84  9.39 

4.89  8.28  5-84  9.59 

4.90  8.30  5.85  9-39 
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h 

0 

3 

4 

8 

12.' 

5 

22 

.5 

32. 

>5 

in  km 

J? 

_s 

_P 

_s 

_P 

J3 

_P 

J3 

_S 

_P 

_S 

_? 

_s 

19.1 

4.00 

7.96 

4.96 

8.03 

4.32 

6.59 

4.47 

7.25 

4.83 

8.62 

4.91 

8.32 

5.86 

9.39 

19.2 

.92 

8.00 

4.98 

3.09 

4.34 

6.59 

4.49 

7.27 

4.85 

8.65 

4.92 

8.33 

5.86 

9.39 

190 

4.95 

8.04 

3.0C 

8.13 

4.36 

6.60 

4.51 

7.29 

4.86 

8.68 

4.93 

8.35 

5.87 

9.40 

1Q.4 

4.97 

8.08 

5.03 

8.17 

4.38 

6,60 

4.53 

7.3I 

4.88 

8.71 

.4.94 

8.37 

5.87 

9.  to 

19.5 

5.00 

8.22 

5.06 

8.22 

4.40 

6.6O 

4.55 

7.33 

4.30 

8.74 

4.96 

8.38 

5.88 

9. to 

19.6 

5.03 

8.17 

5.08 

8.26 

4.42 

6,6l 

4.57 

705 

4.92 

8.78 

4.97 

8.4o 

5.89 

9. to 

19.7 

5.05 

8.21 

5.10 

5.30 

4.44 

6.61 

4 . 59 

7.37 

4.94 

8.8l 

4.98 

8.42 

5.89 

9. to 

19.8 

3.08 

8.25 

5.13 

8.34 

4.46 

6.61 

4.6l 

7.39 

4.95 

8.84 

4.99 

$.44 

5*90 

9.41 

19.9 

5.10 

8  2Q 

5.16 

8.38 

4.48 

6.6l 

4.6^ 

7.41 

4.97 

8.87 

5.00 

8.45 

5.90 

9.41 

20.0 

3.13 

8.33 

5.1b 

8.43 

4.50 

6.61 

4.65 

7.43 

4.99 

8.90 

5.01 

8.47 

5.91 

9.4i 

With  these  calculations  completed  the  task  of  locating  epicenters  was  begun. 

The  data  as  presented  in  Table  7  represents  an  initial  correlation  of  arrivals  et  the 
various  stations  based  on  time.  A  number  of  these  eyeball  correlations  failed  to  be 
consistent  when  the  data  were  used  to  locate  an  epicenter. 

Two  graphical  methods  were  used  in  locating  epicenters  at  this  point.  The 
first  method  was  based  on  A  (P-H)  arc  distance  values  for  each  station.  A  trial 
II  (hypoeentral  time)  value  was  computed  with  the  equation: 

P-H  =  1.37  (S-P)  (5) 

Using  this  H  value,  P-H  values  were  calculated  for  each  station  and  the  correspon¬ 
ding  distances,  ^  (P-H)  were  t  ken  from  the  travel  time  tables.  These  distance- 
values  were  then  plotted  on  a  map  of  the  area.  Intersection  of  the  distance  arcs 
at  a  point  determined  the  epicenter.  If  the  distance  arcs  failed  to  intersect' 
at  a  point  the  H  value  was  adjusted  until  such  an  intersection  was  achieved. 

Should  adjustment  of  H  fail  to  effect  an  intersection,  focal  depth  was  adjusted. 

The  second  method  used  was  the  hyperbolic  ^locus  method.  The  time  difference 
between  a  rivals  at  any  pair  of  stations  will  determine  a  hyperbolic  l-'cm.  by 
arbitrarily  choosing  hypocentral  times.  For  a  given  H,  distance  arcs  front  each 
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station  of  a  chosen  pair  intersect  at  a  point  which  is  a  point  on  the 
hyperbolic  locus.  Additional  points  on  the  locus  are  determined  by 
successively  chosen  values  of  H.  A  second  pair  of  stations  gives  a 
second  hyperbolic  locus  and  similarly  for  a  third  pair.  The  hyperbolic 
loci  should  intersect  at  a  point  which  is  the  epicenter.  Both  the 
hyperbolic  method  and  the  arc -distance  methods  require  a  minimum  of 
three  stations  to  determine  an  epicenter.  Table  7  lists  216  separate 
quakes  but  only  sixty-five  of  these  are  recorded  at  three  or  more 
stations.  Of  these  sixty-five  quakes,  epicenters  were  located  for  forty- 
two.  These  quakes  and  their  epicenters  are  listed  in  Table  15  and  shown 
in  Figure  13  • 


TABLE  15  EPICENTERS  KXLAUEA 


Quako  Max.  Ampl. 


Number 

Epicenter  Station 

Distance 

in  mm. 

- 

19°  24.2'  N 

U 

5.8 

1.10 

155°  19.4'  W 

0 

4.6 

•  52 

H  =  21-59-18.4 
h  =  0  km 

M.L# 

12.4 

.05 

5 

19°  24.4'  N 

U 

3-3 

.15 

155°  19-2'  W 

0 

4.5 

.22 

H  =  21-45-09-7 
h  =  0  tan 

M.L. 

12.3 

.03 

7 

19°  24.8'  N 

0 

6.8 

.10 

155°  20.5'  W 

D 

10.0 

•03 

H  =  22-09-19.1 
h  =  0  km 

M.L. 

10.4 

.20 

l6 

19°  24.2'  N 

U 

5.1 

•  30 

155°  20.2'  W 

0 

5-8 

.40 

H  =  25-50-54.4 
h  =  0  km 

M.L. 

n.7 

.07 
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155*15’ 


>°20 ' 

rl  I  1  1  19°50» 


155  20' 


155*15 ’ 


FIGURE  13  EPICENTERS,  KIUUEA 
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Quake 


Number 

Epicenter  Station 

Distance 

22 

19°  25-0'  N  U 

3-3 

155°  19-4'  W  D 

11.1 

H  =  00-19-25-4  M.L. 
h  =  0  km 

11.4 

25 

19°  24.7'  N  U 

3-4 

155°  19-4'  W  d 

10.8 

H  =  00-28.05.5  M.L. 
h  =  0  km 

11.6 

28 

19°  22.9'  N  U 

5.2 

155°  18.9'  W  0 

3.5 

H  =  00-52-51.2  M.L. 
h  =  0  km 

15.O 

29 

19°  24.0'  N  U 

6.0 

155°  14.2'  W  0 

5.4 

H  =  00-55-45.5  M.L. 
h  =  5*0  km 

18.4 

50 

19°  27-7*  N  U 

6.2 

155*15.0'  W  0 

8.5 

H  =  00-56-07.2  M.L. 
h  =  0  km 

15.1 

34 

19°  25.9*  N  U 

2.8 

155°  18.0 '  w  0 

2.2 

H  =  00-42-05.6  M.L. 
b  =  5-0  km 

14.3 

4l 

19°  25.1'  N  U 

2.2 

155°  18. 8*  W  D 

12.0 

H  =  00-48-50.5  M.  L. 
h  =  0  km 

11.6 

43 

19°  24.9'  N  U 

3-0 

155°  19.2'  W  D 

11.2 

H  =  00-51-52.5  M.L. 
h  =  0  km 

11.7 

46 

19°  24.6'  N  U 

3.3 

155°  19-3'  W  0 

4.6 

H  =  00-55-55.1  M.L. 
h  =  0  km 

12.0 

Max.  Ampl. 
in  mm. 


.20 

.07 

.11 


4.4o 

.19 

.15 


1.40 

•45 

.17 


•  50 
,45 
•05 


1.80 

•55 

.26 


l.6o 

.45 

.12 


9.00 

.22 

•50 


6.00 

.10 

.15 


6.20 

.20 

.10 
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Quake 


Number 

Epicenter 

Station 

Distance 

59 

19°  24.5*  N 

U 

3.1 

155°  19.0*  W 

D 

10.8 

H  =  01-14-37-7 
h  =  0  km 

M.L. 

12.4 

67 

19°  25-3'  N 

U 

1.6 

155°  18. 6*  W 

D 

12.2 

H  =  01-41-03-8 
h  =  0  km 

M.L. 

11.5 

70 

19°  24.0*  N 

U 

3-3 

155°  18. 7'  W 

0 

3-3 

H  =  01-43-32.0 
h  -  0  km 

M.L. 

13.4 

74 

19°  23*9*  N 

U 

3-0 

155°  17.8*  W 

D 

11.7 

H  =  01  49-55.6 
h  =  4.0  km 

M.L. 

l4.6 

81 

19°  25-0*  N 

U 

1-9 

155°  18.5'  W 

0 

4.0 

H  =  02-05-21.7 

D 

12.0 

h  =  4.0  km 

M.L. 

12.0 

84 

19°  23-71  N 

U 

3.2 

155°  16.5'  W 

0 

1.2 

H  =  02-12-58.4 

D 

13.7 

h  =  4.0  km 

M.L. 

16.2 

89 

19°  2k. y  N 

U 

2.3 

155°  18. 5'  W 

0 

3.4 

H=  02-26-42.8 
h  =  0  km 

M.L. 

12.9 

91 

19°  25.8’  N 

U 

3.1 

155°  15- 9'  W 

D 

16.6 

H-02-33-45.2 
h  =  8.0  km 

M.L. 

15.0 

95 

19°  24. y  N 

U 

4.0 

155°  19.5*  W 

0 

5.0 

H  =  02-37-58.1 

D 

10.3 

h  =  0  km 

M.L. 

12.1 

Max.  Ampl. 
in  mm. 

7-50 

•35 

.4o 


3.00 

.10 

.32 


.70 

.20 

.05 


2.00 

.05 


5.70 

1.00 

.36 

.82 

1.10 

1.20 

•  9'. 

.14 

1.80 

.15 

.20 


16.00 
•  30 
.10 


3.00 

•35 

.06 

.16 
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Quake  Max.  Ampl. 


Number 

Epicenter 

Station 

Distance 

in  mm. 

97 

19°  24.4’  N 

U 

2.2 

1.50 

155°  18.1'  W 

:  2.8 

.... 

H  =  02-41-24.2 

I) 

:n.9 

.02 

h  =  3-0  km 

M.L. 

13.5 

.04 

101 

19°  25.3'  N 

0 

4.7 

•  70 

155°  18.8'  W 

N.P. 

3-2 

.20 

H  =  02-48-16.2 

D 

12.1 

.06 

h  =  0  km 

M.  L. 

11.6 

.06 

103 

19°  24.1*  N 

0 

2.7 

.40 

155°  18.3'  W 

N.P. 

2.7 

.80 

H  =  02-57-16.0 
h  =  0  km 

D 

11.3 

•  50 

104 

19°  26.8*  N 

U 

6.9 

I.60 

155°  13-9'  W 

0 

8.2 

.80 

H  =  03-00-20.8 
h  =  3*0  4m 

M.L. 

17.5 

.05 

n6 

19°  25.0'  N 

U 

4.0 

6.50 

155°  15.4'W 

0 

4.4 

.25 

H  =  03-23-23  .'6 

D 

16.6 

•  15 

h  *  0  km 

M.L. 

16.5 

.15 

117 

19°  25.4*  N 

U 

6.7 

15.80 

155°  13-7’  W 

D 

19.2 

•  35 

H  =  03-24-12.4 
h  -  8.0  km 

M.L. 

18.7 

.30 

120 

19°  25-3'  N 

U 

2.6 

12.00 

155°  19-0*  W 

D 

11.9 

.45 

H  =  03-28-30.1 
h  =  0  km 

M.L. 

ll.l 

.60 

125 

19°  23.9’  N 

U 

2.7 

3.00 

155°  17-5'  W 

D 

12.1 

.80 

H  =  03-40-51.2 
h  =  8.0  km. 

M.L. 

14.9 

.10 

135 

19°  23.4'  N 

U 

4.4 

3.50 

155°  18.9'  W 

D 

9.6 

.76 

H  =  04-02-47.6 
h  =  3*0  km 

M.L. 

14.1 

.16 
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Quake  Max.  Ampl. 


Number 

Epicenter 

Station 

Distance 

in  mm. 

144 

19°  24.8*  N 

U 

2.9 

3-50 

155°  19.1'  W 

D 

11.1 

-13 

H=04-l5-24.3 
h  =  0  km 

M.L. 

11.9 

.16 

147 

19°  25.6'  N 

U 

6.4 

6.10 

155°  13-9'  W 

N.P. 

5.6 

8.00 

H=o4-23-19.o 

D 

19.0 

.11 

h  =  4.0  km 

M.L. 

18.1 

.25 

150 

19°  23-6'  N 

U 

3-3 

2.30 

155°  17.5*  W 

D 

11.8 

-31 

H=o4-26-44.2 
h  =  5*0  km 

M.L. 

15-3 

.20 

158 

19°  25.2'  N 

U 

1-7 

36.00 

155°  18.5*  W 

0 

4.2 

— 

H=04-39-15-6 

D 

12.4 

.80 

h  =  3*0  km 

M.L. 

12.0 

1.60 

159 

19°  24.8'  N 

U 

2.1 

l4.oo 

155°  18. 6*  W 

0 

3-8 

-55 

H=o4-43-01.2 

D 

11-9 

.09 

h  =  3.0  km 

M.L. 

12.3 

-77 

168 

19°  25.8'  N 

U 

4.0 

6.50 

155°  15-3'  W 

N.P. 

3-2 

13.00 

H=04-52-31.1 

D 

17.4 

.21 

h  =  5*0  km 

M.L. 

16.0 

.12 

lf4 

19°  24.9'  N 

U 

3.6 

2.70 

153°  19.6'  W 

D 

10.9 

.40 

H=05-05-02.2 
h  =  3*0  km 

M.L. 

11.3 

.18 

181 

19°  25-0'  N 

U 

3-9 

-70 

155°  19-8 1  W 

D 

10.9 

-05 

H=05 -18-35-1 
h  =  0  km 

M.L. 

10.9 

-09 

182 

19°  25-4'  N 

U 

4.0 

1.70 

155°  15-2'  W 

D 

16.9 

-07 

H=05-22-35-8 
h  =  5-0  km 

M.L. 

16.3 

.04 
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Quake 

Number 

Epicenter 

Station 

Distance 

Max.  Ampl. 
in  mm. 

196 

19°  24.4*  N 

U 

1-9 

2.80 

155°  17.3'  W 

D 

15.0 

•  91 

H=05-47-05.9 
h  =  8.0  km 

M.L. 

14.6 

.20 

206 

19°  25.7'  N 

U 

5.0 

8.00 

155°  19.4*  W 

D 

12.4 

.90 

H=06 -09-51.0 
h  =  0  km 

M.L. 

10.4 

.18 

A  number  of  difficulties  arose  in  the  epicenter  program.  There  were  a 
number  of  instances  where  the  data  were  incompatible  and  no  epicenter  could 
be  located.  In  these  instances  arrival  times  at  the  various  stations  were 
not  from  the  same  quake  but  were  arrivals  from  different  earthquakes  that 
occurred  almost  simultaneously.  If  arrivals  from  both  events  reached  a  par¬ 
ticular  station,  one  of  the  arrivals  was  overridden  and  lost.  Arrivals  fre¬ 
quently  failed  to  reach  a  particular  station  because  the  energy  was  atten¬ 
uated  and  absorbed  due  to  its  high  frequency  character  and  low  energy  at  the 
source. 

This  latter  condition  was  very  noticeable  at  the  North  Pit  Station  where 
the  number  of  quakes  recorded  was  much  larger  than  at  any  other  station.  In 
spite  of  this,  there  were  very  few  instances  where  the  North  Pit  data  were 
compatible  with  that  of  tjie. other  stations.  It  also  seems  likely  that  there 
may  be  a  local  velocity  anomaly  at  North  Pit  which  effects  the  arrivals. 

In  a  number  of  cases  it  was  not  possible  to  achieve  a  unique  epicenter. 
By  varying  the  focal  depth  an  epicenter  could  be  found  at  any  one  of  several 
depths  varying  from  the  surface  on  down.  As  the  focal  depth  changed  the 
epicenter  shifted.  This  naturally  gave  rise  to  the  question,  which  epicenter 
to  use?  If  there  was  nothing  in  the  character  of  the  seismogram  to  indicate 
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focal  depth,  then  the  surface  focus  was  accepted  as  the  proper  one  on  the 
"basis  that  these  seismic  events  are  associated  with  the  volcanic  action 
which  was  within  a  matter  of  hours  of  erupting.  Hence,  the  activity  was 
certainly  very  near  to  the  surface. 

7-5  THE  COMPUTES  PROGRAM 

While  the  graphical  location  of  epicenters  was  going  on,  a  computer 
program  was  being  developed  to  process  the  data.  The  program  is  based  on 
the  assumption  of  rectilinear  paths  for  the  seismic  waves  since  the  dis¬ 
tances  are  very  small.  The  basic  equation  of  the  program  is  the  distance 
formula: 

(x  -  x^2  +  (y  -  yt)2  +  (z  -  2i)2  =  V2  (Pt  -  H)2,  i=l,2,3A  (7) 

Where  x^^  -  latitude  of  the  station 
yi  *  longitude  of  the  station 
=  elevation  of  the  station 

V  =  velocity  of  the  longitudinal  wave 
P^  =  observed  arrival  time  of  P  wave 

H  =  hypocentral  time 

For  local  quakes  a  suitable  origin  point  is  selected  to  encompass  the 
network  of  observing  stations  and  the  possible  epicentral  region.  For  the 
earthquakes  under  investigation  in  Hawaii,  the  origin  point  selected  was 
19°  20'  N  and  155°  10s  W.  Latitude  and  longitude  for  each  station  lo¬ 
cation  were  then  converted  to  distance  in  kilometers  from  the  origin  using 
the  conversion  factors  for  one  minute  of  arc  of  meridian  and  one  minute  of 
arc  of  parallel  at  this  latitude.  Station  elevation,  given  in  meters  above 
sea  level,  is  converted  to  kilometers. 
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Equation  7  involves  four  unknowns  for  which  ve  wish  to  solve  x,  y,  z, 
and  H.  Pour  equations  are  required  so  that  observations  at  four  stations 
are  necessary  for  a  solution.  In  order  to  program  the  solution  it  was 
necessary  to  express  the  equations  in  linear  fora.  A  mathematical  reduc¬ 
tion  of  the  four  Equations  7  by  eliminating  the  variable  H  results  in  a 
linear  system  of  three  equations  of  the  form: 

A^x  +  B^y  +  C^z  +  =  0  i  =  1,2,3  (8) 

Where  A±  =  2  [(P4  -  P2)(x4  -  x3)  -  (P4  -  P3)(x4  -  x2)] 

\  =  2  [ (P4  -  P2)(y4  -  y3)  -  (P4  -  P3)(y4  -  y2)] 

=  2  [(P4  -  P2)(z^  -  z3)  -  (P4  -  P3)(z4  -  z2)] 

Dx  =  (P4  -  P2)(x32  +  y32  +  z32) 

(p4  -  r3)(x22  +  y22  +  z22) 

+  (PjB  -  P2)(x42  +  y42  +  z42)  ' 

+  v2  (p4  -  P3)(P4  -  P2)(P3  -  p2) 

Similar  values  of  the 
coefficients  occur  for 
A2  B2  C2  I>2  and  A3  B3  C3  D3 

Unfortunately  the  system  7*52  is  not  an  independent  system.  Any  pair 
of  Equations  8  is  independent  but  the  third  equation  is  always  dependent. 
Geometrically  each  Equation  8  represents  a  plane  and  if  the  system  were  in¬ 
dependent  the  solution  would  be  a  point  which  represents  the  epicenter. 
Since  the  system  is  dependent  the  solution  will  be  a  line  solution. 

When  solutions  of  Equations  8  were  attempted  the  results  were  erratic. 
This  was  disturbing  but  since  the  equations  represented  three  planes  deter¬ 
mined  by  the  station  coordinates  and  the  observed  arrival  times,  it  T*as  at 
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first  felt  that  the  result  might  be  legitimate.  The  coelf icients  are  de¬ 
termined  by  the  observational  data  which  is  limited  in  its  accuracy  and  the 
resulting  numerical  values  for  the  coefficients  did  not  reveal  the  depen¬ 
dent  neture  of  the  equations.  A  reexaminabion  of  the  theoretical  solution 
of  Equations  8  revealed  the  dependent  nature  of  the  system. 

This  meant  that  our  solution  set  of  eguations  consisted  of  two  equa¬ 
tions  in  three  unknowns,  so  a  method  had  to  be  devised  to  handle  this.  The 
three  unknowns  in  Equations  8  are  x,  y,  and  z.  The  method  consisted  of  the 
following . steps : 

1.  Assume  a  value  for  z. 

2.  Substitute  z  in  Equations  8  and  compute 
the  values  of"x  and  y. 

3.  Substitute  the  assumed  val’  -.  of  z  and  the 
computed  values  of  x  and  y  in  the  four 
Equations  7  and  compute  a  value  of  H  for 
each  equation. 

If  the  four  values  of  H  agree  then  the  set  of  values  x,  y,  z,  and  H  are  the 
solution.  The  usual  circumstance  is  that  the  H  values  will  not  agree.  If 
they  do  not  agree  then  a  new  value  is  assumed  for  z  and  the  process  repeated . 
In  practice,  values  of  z  were  chosen  for  every  0.5  kilometers .  As  the  z 
values  are  varied,  the  H  values  for  each  station  will  converge  to  a  common 
value  which  represents  the  solution.  As  this  convergence  occurs,  the  in¬ 
dividual  values  of  H  for  a  given  equation  will  rise  through  a  maximum  value 
and  then  fall  off.  This  maximum  value  is  the  convergence  value.  The  full 
import  of  this  rise  to  a  maximum  is  not  clear  at  this  time. 

This  method  worked  very  well  and  yielded  epicentral  solutions  which 
gave  an  independent  check  on  the  graphical  solutions  previously  achieved. 
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It  should  be  pointed  out  that  in  the  computer  solution,  values  for  2  the 
focal  depth,  can  legitimately  be  positive  since  the  elevation  of  the  Kilauea 
Crater  is  about  one  kilometer. 

In  the  graphical  solutions  when  difficulty  was  experienced  in  obtain¬ 
ing  an  epicentral  fix  due  to  an  incorrect  P  arrival  at  one  of  the  stations, 
it  was  not  always  clear  which  station  was  in  error.  However,  in  the  com¬ 
puter  program  if  there  is  an  incorrect  value,  this  will  stand  out  '^ry  clear 
ly  in  the  print  out  result,  since  the  legitimate  values  will  produce  conver¬ 
gence  to  a  common  value  for  H  but  the  incorrect  value  will  not. 

This  result  led  to  a  further  adaptation  of  the  computer  program  to 
quakes  which  were  recorded  at  only  three  stations.  When  an  incorrect  P  ar¬ 
rival  time  value  occurred  which  failed  to  produce  a  convergence  for  H,  this 
incorrect  arrival  was  adjusted  until  its  H  value  did  converge.  Obviously, 
this  adjusted  arrival  time  is  fictitious  and  does  not  exist.  For  the  quakes 
for  which  only  three  stations  had  P  arrivals,  a  fictitious  P  arrival  was 
fabricated  for  the  fourth  station.  This  set  of  data  with  the  one  spx’rious 
value  was  put  into  the  computer,  and  this  spurious  value  was  adjusted  until 
the  computer  produced  a  convergent  solution.  This  procedure  is  not  as 
arbitrary  as  it  might  seem.  Basically  it  amounts  to  asking  the  computer  to 
find  a  fourth  P  arrival  that  will  fit  the  other  three  observed  arrivals, 
given  the  crustal  structure,  velocity  and  coordinates  of  the  stations. 
Another  way  of  looking  at  it  is  that  normally  the  computer  is  supplied  with 
sixteen  pieces  of  information  and  asked  to  compute  four  others.  In  the  new 
scheme  the  computer  is  supplied  with  only  fifteen  pieces  of  information  and 
asked  to  compute  five  others. 
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As  a  result  of  this  we  have  been  able  to  check  the  graphical  epicenters 
and  the  computer  epicenters  against  each  other  and  find  that  there  is  excel¬ 
lent  agreement. 

7.6  FIRST  MOTION  PATTERNS 

§y  combining  the  firs  „  motion  data  with  the  epicenter  data  the  first 
motion  pattern  for  each  quake  was  obtained.  The  data  showed  a  fairly  high 
consistency,  with  three  patterns  resulting.  When  the  compression  -  dilata¬ 
tion  data  were  plotted  at  each  station,  quadiantal  lines  were  drawn  through 
the  epicenter  to  separate  the  areas  of  dilatation  from  the  areas  of  com¬ 
pression.  In  scxne  cases  the  position  of  these  lines  could  be  varied  by  as 
much  as  ten  to  fifteen  degrees  of  rotation  about  the  epicenter  but  this  did 
not  fundamentally  change  the  pattern.  In  other  cases  the  pattern  was  rigid¬ 
ly  fixed. 

One  of  the  quadrantal  lines  shows  a  trend  Northeast  -  Southwest,  and 
the  other  at  right  angles  to  this  shows  a  Northwest  -  Southeast  trend.  The 
Northeast  -  Southwest  direction  coincides  with  the  direction  of  the  South¬ 
west  Rift  Zone  and  the  coincidence  of  these  directions  probably  represents  a 
relationship  between  this  tectonic  feature  and  the  seismic  motion. 

If  we  assuuia  that  this  Northeast  -  Southwest  direction  represents  the 
line  along  which  the  seismic  motion  occurs  then  in  fifteen  cases  the  North¬ 
west  side  of  the  block  shifted  to  the  Southeast.  In  fourteen  cases  the 
Northwest  side  shifted  to  th'  Northeast.  These  represent  two  of  the  three 
main  trends  observed.  They  represent  opposite  direction  of  motion  and  seem 
almost  equally  likely  to  occur. 
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As  the  epicenters  shift  from  the  area  Southwest  of  Kilauea  Caldera  to 
an  area  Northeast  of  it,  the  epicenters  tend  to  get  deeper  and  the  trend  of 
the  quadrantal  line  shifts  from  a  Northeast  -  Southwest  direction  to  a  dir¬ 
ection  more  East  -  West.  Here  again,  a  general  correlation  can  be  made  bet¬ 
ween  this  direction  and  that  of  the  East  Rift. 

Die  coincidence  of  these  quadrantal  lines  with  the  Southwest  Rift  and 
with  the  East  Rift  seems  more  than  accidental  and  augurs  for  a  mechanistic 
connection.  However,  if  the  seismic  motion  occurs  along  the  other  quadran¬ 
tal  line  which  trends  Northwest  -  Southeast,  motion  along  the  Southwest 
part  o"  the  block  would  be  Northwesu  or  Southeast. 

In  seven  cases  horizontal  motion  could  not  explain  the  distribution  of 
motion.  For  these  a  vertical  motion  was  assumed  with  two  blocks  separated 
by  a  Northwest  -  Southeast  trending  line.  The  Northeast  block  was  char¬ 
acterized  by  upward  motion  and  the  Southwest  block  by  downward  motion.  The 
three  patterns  of  motion  are  shown  in  Figure  14,  15,  and  16. 


7. 7  AMPLITUDE  -  DISTANCE  REIATIONSHIPS 

Maximum  amplitude  values  at  each  station  for  each  earthquake  were  com¬ 
bined  with  the  corresponding  distances  computed  in  the  epicenter  program. 

Die  medial  test  was  used  to  determine  whether  or  not  a  significant  relation¬ 
ship  existed  between  amplitude  and  distance.  It  will  be  recalled  that  in 
the  case  of  Vesuvian  quakes  no  significant  relationship  was  found.  The 
medial  test  is  shown  in  Figure  17« 

On  the  strength  of  the  medial  test  indicating  a  relationship  between 
amplitude  and  distance  an  investigation  into  magnitude  for  these  Hawaiian 
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155°  20' 


FIGURE  14  EARTH  MOTION  KHAUEA,  PATTERN  I 
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155°20' 


155 *15' 


FIGURE  15  EARTH  MOTION  KELAUEA,  PATTERN  II 
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155 0  20 f  1^5°  IV 


155*  20'  155°  1 

FIGURE  1 6  EARTH  MOTION  KXLAUEA,  PATTERN  III 
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volcanic  quakes  was  begun.  Richter's  relationship  for  shallow  earthquakes 
in  California  between  magnitude  M,  maximum  trace  amplitude  A,  and  distance 
is 

M  =  log  A  +  3  log  A  -  3.37  (9) 

When  data  for  individual  quakes  were  plotted  as  log  A  versus  log  A  , 
the  points  tended  to  scatter,  slope  values  were  difficult  to  determine  and 
showed  a  wide  variation.  At  this  time  it  appears  that  the  use  of  data  from 
a  single  quake  is  not  reliable.  Conclusions  based  on  a  single  event  should 
be  avoided  until  individual  station  corrections  have  been  computed. 

The  data  for  all  the  quakes  of  known  epicenters  were  grouped  together 
on  a  single  graph  of  log  A  versus  log  A.  The  result  was  a  fairly  homo¬ 
geneous  picture.  A  center  line  was  eyeballed  through  the  data  and  the  equa¬ 
tion  of  this  regression  line  was  computed  to  be 

log  A  +  2.3  log  A  =  1.64  (10) 

Expressing  magnitude  M  as  a  function  of  A  and  A  gives  the  relationship 

M  =  log  A  +  2.3  log  A  -  1.64  (11) 

Applying  Richter's  definition  of  a  zero  magnitude  earthquake,  namely  M  =  0, 

A  =  10  3mm,  A  =  100  Ion  to  Equation  11  indicates  that  the  quantities  deter¬ 
mined  empirically  from  the  observed  data  are  in  reasonable  agreement  with 
the  magnitude  concept. 

For  an  earthquake  of  given  magnitude,  Richter's  equation  indicates 
that  the  amplitude  varies  inversely  as  the  cube  of  the  distance  in  Southern 
California.  For  Hawaiian  quakes,  Equation  11  indicates  that  for  a  shock  of 
given  magnitude,  the  amplitude  varies  inversely  as  the  2.3  power  of  the 


74 


distance.  Thus,  amplitude  seems  to  attenuate  less  with  distance  in  Hawaii 
than  in  California.  This  may  he  correlated  with  the  folded  and  faulted 
nature  of  the  rocks  in  the  California  region.  It  is,  however,  somewhat  sur¬ 
prising  since  high  attenuation  does  seem  to  occur  in  the  Kilauea  Summit 
Region,  particularly  in  the  area  of  the  North  Pit  Station  which  lies  in  the 
caldera.  This  fact  probably  accounts  for  a  large  part  of  the  attenuation 
since  the  rocks  in  the  caldera  region  are  probably  weak.  The  high  frequency 
character  of  the  volcanic  -  seismic  events  probably  also  accounts  in  part 
for  the  attenuation.  These  factors  seem,  therefore,  to  have  a  pronounced 
local  effect  which  drops  off  with  the  distance. 

Using  Equation  11,  curves  for  magnitude  M  =  +1,  M  =  0,  and  M  =  -1  were 
computed  and  plotted  on  the  graph  of  the  log  A  versus  log  £  .  Practically 
all  the  quakes  fell  between  these  magnitude  values  which  indicates  the  very 
low  energy  level  of  such  shocks.  The  graph  of  these  data  is  shown  in 
Figure  18. 
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0.1  1.0  10 
Distance  (km) 


FIGURE  18  LOG  A  VS  LOG  b,  KELAUEA 


76 


8 


JAPANESE  VOLCANOES 


The  Japanese  Islands  possess  many  volcanoes  and  these  have  become  the 
object  of  intensive  research.  Many  of  the  volcanoes  are  explosive  in  char¬ 
acter  and  excellent  observatories  covering  practically  all  aspects  of  geo¬ 
physical  science  applicable  to  volcanoes  have  been  set  up  to  investigate 
the  volcanic  phenomena. 

One  such  observatory  is  located  at  Mount  Asama.  Through  the  courtesy 
and  kindness  of  Professor  Ryutaro-Takahasi,  Director  of  The  Earthquake 
Research  Institute,  The  University  of  Tokyo,  preliminary  arrangements  were 
made  to  visit  the  Institute  and  examine  the  seismic  data  from  Mount  Asama. 
Final  arrangements  were  completed  through  the  kindness  of  the  new  Director, 
Professor  Hirosi  Kawasumi.  Professor  Takeshi  Minakami,  Director  of  The 
Mount  Asama  Observatory  and  Professor  of  Physical  Volcanology  and  his  As¬ 
sistant,  Dr.  Kiyoo  Mogi,  rendered  invaluable  assistance  in  making  the  seis¬ 
mic  records  available  for  analysis. 

8.1  MOUNT  ASAMA  VOLCANO  OBSERVATORY 

Mount  Asama  is  located  approximately  one  hundred  kilometers  west  and 
slightly  north  of  Tokyo.  The  Asama  Volcano  Observatory  is  located  on  the 
eastern  flank  of  the  volcano  approximately  four  kilometers  from  the  summit 
crater.  Asama  is  one  of  the  most  active  of  the  Japanese  volcanoes  and 
eruptions  are  frequently  explosive.  Such  an  explosive  eruption  occurred  on 
August  18,  1961.  The  seismic  activity  associated  with  this  event  is  the 
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material  investigated  here. 

At  the  time  of  the  1961  eruption  a  seismic  network  of  seven  stations 
was  operating.  Seismic  signals  are  telemetered  from  the  instruments  to  the 
Asama  Observatory  where  they  are  recorded  on  smoked  paper  which  gives  a 
very  fine  definition  of  the  high  frequency  motion.  The  seismographs  are 
elect romagnetic  and  are  oriented  as  vertical  or  horizontal  radial  compo¬ 
nents.  The  station  designations  and  seismographs  located  p.t  each  are  given 
in  Table  l6. 

TABLE  16  MOUNT  ASAMA  SEISMOGRAPH  STATIONS 


Station 


Coordinates 


North 

East 

Z  (km) 

A 

HR 

36°  23-9' 

138°  32.2* 

2.2 

B 

HR,Z 

36°  24.1' 

138°  32.9 1 

1.8 

C 

HR,Z 

3 6°  24.0* 

138°  34.0* 

1.4 

D 

HR 

36°  22.6' 

138°  32.2* 

1.5 

E 

HR 

36°  23.3' 

138°  33.3' 

1.6 

F 

HR 

3 6°  25.6' 

138°  32-7’ 

1.5 

G 

Z 

3 6°  23.6' 

138°  29.8’ 

2.0 

HR  =  Horizontal  Radial  component 
Z  =  vertical  component 


In  the  investigations  at  Mount  Asama,  the  Japanese  scientists  are 
primarily  interested  in  the  question  of  prediction  of  volcanic  eruptions. 
Professor  Minakami  has  developed  a  statistical  method  for  such  prediction 
based  on  observation  of  the  number  of  seismic  events  occurring.  The  ques¬ 
tion  of  hypocentral  depth  is  important  and  is  determined  from  the  char¬ 
acter  of  the  seismogram.  Since  the  seismic  events  eire  considered  as  orig¬ 
inating  in  the  volcano  structure,  location  of  epicenters  is  not  carried  out 
and  seismic  velocities  and  local  structure  apparently  have  not  beer,  invest¬ 
igated.  Two  types  of  quakes,  A  and  B,  are  distinguished  by  Minakami. 


A-type  quakes  are  deep  with  the  hypocenters  from  one  to  ten  kilometers .  B- 
type  quakes  are  shallow  with  the  hypocenters  from  zero  to  one  kilometers. 
Earthquakes  of  A-type  are  similar  to  tectonic  quakes  and  show  distinct  P  and 
S  phases.  B-type  quakes  frequently  begin  with  large  motions  and  show  char¬ 
acteristically  large  surface  waves. 

Minakami  refers  to  explosion  earthquakes  as  resulting  from  individual 
explosions  originating  in  the  volcano.  The  hypocenter  of  such  an  explosion 
is  situated  not  much  deeper  than  the  active  crater  floor.  Initial  motion 
is  a  push  in  all  directions. 

8.2  SEISMOGRAPH  READINGS 

Prior  to  the  August  l8th  eruption,  seismic  activity  began  tc  increase. 
Volcanic  explosions  occurred  intermittently  with  increasing  intensity  and 
frequency,  culminating  in  the  paroxysm  at  14  4l  5Cr  local  time  August  18, 
1961.  At  about  1000  hours,  seismic  recording  was  switched  off  because  of 
the  interference  coming  from  electrical  phenomena  associated  with  the  erup¬ 
tion.  Recording  was  resumed  at  approximately  1500  hours  after  the  parox- 
ysmic  eruption. 

The  following  seismic  readings  are  given  only  in  seconds,  with  hours 
and  minutes  omitted.  Direction  of  motion  is  indicated  as  dilatation  or  com¬ 
pression  for  all  components.  For  the  vertical  component  this  has  the  usual 
meaning.  For  the  horizontal  radial  components,  a  push  is  designated  as  a 
compression  and  a  pull  as  a  dilatation. 


79 


TAELE  17  SEISMOGRAPH  READINGS 


Station  _ A 

Quake  _ H 


H 


B 


D 

H 


E  F  G 

H  H  Z 


August  11 

1 

42.8 

c 

^3.3 

c 

2 

12-4 

d 

4o .  5 

August  12 

3 

39-5 

c 

c 

4 

15-5 

c 

15-5 

c 

15.2  d 

15.2  c 

15.1  c 

5 

20. 2 

c 

20.6 

c 

20.2  d 

21.1  c 

21.5  c 

21.4  c 

6 

47.4 

d 

46.6 

c 

7 

57.° 

c 

56.8 

d 

56.9  c 

58.3  c 

8 

33-2 

c 

34.3 

d 

34.2  d 

350  d 

35.3  d 

9 

39-5 

c 

40.8 

c 

10 

51.3 

c 

August  13 

li 

20.2 

d 

20.2 

c 

12 

31.8 

d 

32.0 

c 

13 

03.0 

c 

02-5 

c 

14 

32.2 

d 

15 

49.7 

c 

49.8 

d 

l6 

01.0 

d 

01.1 

c 

01.6  d 

01.3  d 

17 

38.0 

d 

16 

58.5 

d 

19 

53.5 

d 

54.0 

c 

20 

10.0 

d 

18.7 

c 

21 

24.5 

d 

24.7 

c 

27.2  d 

28.3  d 

28.5  c 

27.5  d 

24.8  c 

22 

42.2 

d 

42.5 

c 

44.6  d  44.6  c 

45.0  c 

42.7  d 

•4.2  c 

42.2  c 

23 

02.2 

c 

03.I 

d 

24 

15.7 

d 

25 

43.0 

c 

42.9 

d 

45.0  d 

43.4  c 

44.1  c 

26 

01.7 

d 

02.3 

c 

02.7  d 

27 

00.4 

d 

01.3 

d 

01.8  c 

28 

54.1 

d 

55-2 

d 

29 

31.0 

c 

32.0 

c 

33.5  d 

33.0  c 

33.6  d 

34.3  d 

30 

05-8 

d 

31 

10-3 

c 

32 

19-2 

d 

20.2 

d 

35 

45.5 

d 

47.3 

d 

47.3  c 

47.7  d 

3^ 

02.7 

i 

02.6 

d 

03.2  c 

02.1  c 

35 

28.2 

c 

29.5 

c 

3 6 

07.8 

c 

07.9 

d 

08.8  d 

08.9  c 

09.3  d 

09.4  d 

10.1  c 

07.5  c 

37 

29.2 

d 

30.2 

d 

29.2  d 

31.3  c 

30.2  c 

29.5  c 

29-7  c 

36 

17.3 

c 

17.3 

d 

18.4  C 

17-5  d 

17.6  d 

19.1  d 

17.5  c 

39 

55.9 

d 

37.7 

d 

80 


:+  ation 

A 

B 

C 

D 

E 

F 

Quake 

H 

H 

Z 

H  Z 

H 

H 

H 

August  14 

40 

47.8  d 

48.0  d 

41 

10.5  d 

11.6  c 

42 

19. C  c 

21.5  d 

43 

18.4  d 

44 

26.9  c 

45 

37-8  d 

39-5  d 

46 

47.3  c 

49.0  c 

47 

06.8  c 

05.0  c 

07-0  c 

05.7  c 

05.7  c 

06.1  d 

48 

51.3  U 

51.0  d 

49 

09.7  d 

09-7  c 

10.4  d 

August  15 

50 

52.6  c 

52.3  d 

53-2  c 

53-7  0 

52.0  d 

51 

06.0  d 

06.8  I 

06.2  c 

07.0  c 

05.8  d 

52 

14.0  d 

15.0  c 

15.1  C 

15.6  d 

53 

37.3  c 

36.7  d 

54 

40.0  c 

55 

20.4  d 

56 

00.5  d 

00.5  d 

00.2  c 

01.0  e 

August  l6 

57 

44.8  c 

45.2  c 

58 

45.5  d 

45.9  d 

46.2  d 

45.6  c 

46.0  d 

46.7  d 

59 

31.2  c 

32.0  c 

60 

45.5  c 

45.5  c 

45.5  d 

45.0  d 

61 

35-3  d 

35-5  C 

35.1  d 

35.2  d 

62 

26.0  c 

26.2  d 

25-9  c 

25.9  d 

63 

55-2  d 

55-5  d 

64 

56.9  C 

58.5  c 

00.9  c 

57.4  d 

57.4  d 

56.5  d 

65 

52.6  d 

52-7  d 

51-7  c 

52.7  c 

66 

30.2  d 

30.9  c 

34.9  d 

33-7  c 

33-7  d 

67 

57-0  c 

58.5  C 

58.4  c 

68 

15.6  d 

17.6  d 

69 

21.5  d 

70 

58.3  d 

71 

22.7  d 

23.1  c 

72 

52.8  d 

52.8  d 

54.0  d 

52.8  c 

52.6  d 

73 

28.4  d 

29.I  d 

30.0  d 

74 

04.8  c 

05.3  c 

75 

18.1  d 

18.7  c 

76 

09.8  c 

77 

52.7  d 

54.7  d 

78 

37-9  c 

38.2  c 

79 

28.5  C 

28.8  d 

29.1  d 

80 

21.8  d 

21.6  c 

22.0  c 

21.8  e 

21.7  d 

81 

12.4  c 

13.5  d 

82 

18.9  c 

20.0  d 

05.5  c 


00.5  d 


56.4  c 

52.4  d 


22.0  d 
52.2  c 


81 


csi  a 


August  17 


83 

06.2  c 

0 6.6  c 

08.0  c 

84 

48.1  c 

48.5  d 

8s 

44.7  d 

86 

12.1  d 

87 

09.6  d 

88 

22.7  d 

23.2  c 

89 

21.1  c 

21.2  a 

90 

02.7  c 

91 

07.5  c 

07.7  c 

92 

22.6  c 

22.8  c 

93 

02.4  d 

94 

50.6  c 

95 

27.3  c 

27-6  d 

96 

25-5  c 

97 

18.3  d 

17.8  d 

98 

40.6  d 

40.8  c 

99 

46.7  c 

100 

44.8  c 

46.5  d 

101 

05.0  d 

05.3  d 

05.0  c 

06.9  d 

102 

47.2  d 

48.0  d 

47.4  c 

49.3  c 

103 

08.7  c 

12.0  c 

11.2  c 

io4 

49.9  c 

50.8  c 

53.0  d 

105 

21.7  c 

22.8  d 

23.9  d 

106 

36.O  c 

39-3  d 

107 

36.8  c 

38.9  d 

39.6  c 

108 

32.6  c 

109 

58.7  c 

59-5  c 

110 

38.8  c 

37-6  c 

37-9  d 

111 

37-3  d 

37-2  c 

38.1  c 

112 

43.1  d 

113 

17.9  c 

19.4  c 

114 

12.8  d 

12.8  d 

115 

01.0  c 

01.2  c 

116 

22.6  d 

24.0  c 

22.8  c 

117 

22.8  c 

118 

05-5  c 

05.9  c 

119 

13.6  d 

14.0  c 

120 

38.5  d 

37-7  c 

121 

20.1  d 

20.7  c 

122 

24.6  c 

123 

38.9  d 

124 

40.1  d 

42.4  c 

125 

11.3  c 

12.2  c 

126 

47.1  d 

47.4  c 

47.2  d 

47.6  c 

127 

11.2  d 

12.8  d 

128 

21.6  d 

22.5  d 

129 

10.5  d 

11.9  c 

11.0  c 

15.9  c  14.0  d 


07.8  c 


29.0  c 

28.2  d  28.3  d  27.2  c 

19.7  d  20.0  c  20.7  c 

40.7  c 

46.5  c  48.2  c  49.O  c 

05.6  d  05.5  d  06.1  d  04.6  d 

48.7  d 

11.2  d  12.5  d  13.0  c  11.3  c 

51.7  d  51.6  d 

23.5  c  23.6  c  23.5  d 

39. C  d  38-7  c  38-8  d 

39-1  c  38.8  c 

00.0  c  59-1  c 

38.6  d 


13.3  c  13.5  c  12.6  c 
24.5  c  24.1  c  25.0  c 

05.9  c 


20.8  c 


11.7  c  12.0  c 
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Station 

A 

B 

C 

D 

E 

F 

Quake 

H 

H 

_ 

Z 

H 

Z 

H 

H 

H 

150 

25.7 

c 

26.0  d 

26.2  d 

25.7  c 

25.9  c 

26.5  d 

151 

46.7  d 

46.7 

c 

47.0  c 

47.0  d 

47.5  c 

47.6  d 

47.5  c 

47.5  c 

152 

17.8  d 

155 

55-6  d 

154 

55.0  c 

56.5 

c 

55-8  c 

56.6  c 

57.5  c 

56.5  c 

56.6  c 

155 

27-5  d 

156 

47.6  c 

157 

42.6  c 

42.8 

c 

42.5  c 

158 

29.2  c 

50.2 

d 

159 

20.5  C 

i4o 

57.7  C 

l4i 

29.0  C 

142 

55.7 

c 

55.2  d 

54.7  c 

55-5  c 

56.0  c 

145 

19.4 

c 

144 

51.7  C 

145 

56.8  d 

56.9 

c 

57.6  d 

146 

15.6  c 

15.8 

c 

14.8  c 

15.6  c 

147 

16.1  d 

17.1  d 

148 

18.4  c 

18.2 

d 

19.2  c 

19-6  c 

19.1  d 

August  18 

149 

05.2  c 

150 

07.2  c 

151 

12.8  c 

152 

58.4  c 

155 

08.8  d 

09.I 

c 

154 

16.5  d 

16.1 

d 

16.2  d 

l6.6  d 

16.4  d 

16.0  c 

16.2  : 

l6.4  c 

155 

45.2  c 

45.8 

c 

45.1  c 

47.5  d 

156 

54.8  d 

58.1  d 

157 

51-5  c 

51.5 

c 

158 

l6.4  c 

17-4 

c 

17.0  d 

17.6  c 

17.5  c 

159 

25.9  C 

24.5 

d 

24.5  c 

25-7  d 

160 

48.2  c 

49.6 

c 

The  total  number  of  push  and  pull  motions,  labeled  c  and  d  respectively 
in  Table  17,  for  the  horizontal  radial  component  of  motion  were  counted, 
giving  the  results  shown  in  Table  l8. 

TABLE  18  HORIZONTAL  RADIAL  MOTION 


_A_  B_  _C_  _D_  J2_  _F_  _G_ 

Push]  72  64  17  50  24  25  10 

Pull]  68  54  15  12  15  20  11 
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NO 


These  figures  show  that  in  all  cases  the  push  motion  predominates  at  all 
stations  and  for  stations  D  and  E  the  predominance  is  very  large.  However, 
the  significant  point  is  that  the  pull  type  motion  occurs  very  frequently. 
Just  how  many  of  these  quakes  are  volcanic  quakes  and  how  many  are  explo¬ 
sion  quakes  is  unknown,  bdt  on  the  basis  of  the  quakes  examined,  it  is  ev¬ 
ident  that  the  seismic  motion  is  not  exclusively  push  or  pull. 

8.3  EPICENTERS 

Although  no  crustal  model  or  velocities  were  available  for  the  Mount 
Asama  Region,  it  was  decided  to  assume  values  for  these  quantities  and  at¬ 
tempt  to  locate  epicenters.  A  velocity  to  2.0  kilometers  per  second  was 
used  for  the  computation  and  was  assumed  to  apply  to  focal  depths  of  sev¬ 
eral  kilometers.  Epicenters  were  located  both  graphically  and  by  the  com¬ 
puter  method.  There  were  sixty-six  such  cases  in  the  observed  data.  Forty 
of  these  were  chosen  on  which  to  attempt  epicenter  locations.  Twenty-seven 
of  these  cases  failed  and  only  thirteen  gave  satisfactory  epicenters.  These 
epicenters  are  shown  in  Table  19* 

TABLE  19  EPICENTERS  -  MOUWT  ASAMA 


Quake  Epicenter  Station 

k  36°  23-9’  N  A 

138°  30.0*  E  B 

H  =  13.5  D 

h  =  +1.0  km  F 

T  36°  23-2’  N  A 

138°  33.1'  E  B 

H  =  55-9  D 

h  =  +1.5  km  F 
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Quake  Epicenter  Station 


25  56°  23.9’  N  A 

138°  31.6'  E  B 

H  =  42.2  D 

h  =  +1.0  km 

34  36°  23-5'  N  A 

138°  31.4'  E  B 

H  -  1.1  C 

h  =  0  km  D 

50  36°  23.6'  N  B 

138°  30-7'  E  E 

H  =  51*2  F 

h  =  +1.0  km  G 

51  36°  23-3'  N  B 

138°  31.1'  E  E 

H  =  4.5  F 

h  =  +1-5  Ion-  G 

58  36°  23.6'  N  B 

138°  31.5*  E  D 

H  =  44.4  E 

h  =  0  km  F 

65  36°  25.0*  N  A 

138°  30. O'  E  B 

H  =  50.1  E 

h  =  0  km  G 

109  36°  25-0'  N  A 

138°  30.8'  E  B 

H  =  57-4  D 

h  =  0  km  G 

114  36°  24.6'  N  B 

138°  31. O'  E  D 

H  =  11.0  E 

h  =  -0-5  km  F 

ll6  36°  24.3'  N  B 

138°  30. O'  E  D 

H  =  21.8  E 

h  =  0  km  F 


85 


Quake 

Epicenter 

Station 

150 

3 6°  24.4*  N 

B 

138°  30.0'  E 

D 

H  =  230 

E 

h  =  0  km 

F 

142 

36°  23.1'  N 

B 

138°  33.3*  E 

D 

H  =  54.8 

E 

h  =  +0.5  km 

F 

These  locations  were  plotted  individually  with  first  motion  data  for  the 
horizontal  radial  component.  The  results  were  not  very  satisfactory.  In 
five  cases  a  quadrantal  pattern  of  motion  could  be  observed  but  the  pat¬ 
terns  were  not  all  similar.  There  were  four  other  cases  of  push  and  pull 
motions  for  which  no  quadrantal  pattern  could  be  observed.  Finally,  there 
were  four  cases  in  whi^h  all  motion  was  push  or  compressional,  three  of 
which  were  logical  patterns.  In  the  fourth  case,  the  pattern  of  motion 
could  not  be  reconciled  with  the  location  of  the  epicenter. 

The  seismic  stations  at  Mount  Asama  are  all  located  on  the  flanks  of 
the  volcano  at  elevations  below  the  summit.  Hypocenters,  therefore,  may 
be  either  above  or  below  the  elevation  of  the  seismic  station.  Interpre¬ 
tation  of  the  earth's  motion  must  take  into  consideration  the  epicentral 
location,  focal  depth,  and  the  horizontal  and  vertical  components  of  seis¬ 
mograph  motion.  If  the  motion  originates  above  the  seismic  station  and  is 
a  push,  it  will  appear  as  a  push  on  the  horizontal  radial  seismograph  and 
as  a  dilatation  on  the  vertical  seismograph  since  the  vertical  component  is 
downward.  The  same  seismograph  motions  would  result  if  the  motion  orig¬ 
inated  as  a  pull  below  the  station  level.  In  this  latter  case,  without 
knowledge  of  the  epicentral  location  and  focal  depth,  the  motion  could 
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easily  be  misinterpreted  as  a  push  originating  in  or  near  the  summit  crater. 
In  circumstances  such  as  this,  interpretation  of  ground  motion  becomes 
rather  complex.  The  assumptions  made  concerning  the  P  wave  velocity  and  a 
crustal  model,  present  additional  complications  which  may  account  for  some 
of  these  interpretive  difficulties .  An  example  of  the  intex"pretation  of 
seismic  motion  and  earth  motion  is  shown  in  Figure  19* 

8.4  FREQUENCY  OF  OCCURRENCE  OF  QUAKES 

One  of  the  outstanding  characteristics  of  volcanic  earthquakes  is  the 
fact  that  they  occur  in  swarms  preceding  an  eruption.  Numerous  investiga¬ 
tions  of  earthquake  frequency  have  been  made.  In  general,  the  smaller  the 
magnitude  of  the  shock,  the  greater  the  frequency  of  occurrence.  Gutenberg 
and  Richter  in  SEISMICITY  OF  THE  EARTH  show  that  the  frequency  of  occurrence 
N  and  the  magnitude  M  computed  from  surface  wave  amplitude  are  related  by 
the  equation 

log  N  =  a  +  b  (8  -  M)  (12) 

Where  a  and  b  are  constants 

In  Japan,  M.  Ishimoto  and  K.  Iida  have  developed  another  empirical 
equation  relating  the  frequency  of  occurrence  N  and  the  maximum  trace  ampli¬ 
tude  A  of  the  seismogram.  This  equation  is 

NAm  -  c  (13) 

Where  m  and  c  are  constants 

The  magnitude  M  of  Equation  12  is  determined  from  the  maximum  trace 
amplitude.  It  was  shown  by  Suzuki  that  the  b  of  Equation  12  and  the  m  of 
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Equation  13  are  so  related  that 


b  -  m  -  1 


(1*0 


In  practice,  amplitudes  are  read  over  a  range  A  to  A  +  dA  and  the  num¬ 
ber  in  this  group  is  then  tabulated.  In  addition  to  tabulating  the  fre¬ 
quency  or  number  of  occurrences  within  a  range  dA,  the  frequencies  can  be 
summed  up  so  that  a  value  is  obtained  for  the  total  number  of  quakes  produc¬ 
ing  all  amplitudes  from  the  largest  down  to  a  specified  value  of  A. 

Minakami  in  Japan  has  used  this  approach  to  show  that  A-type  and  B-type  vol¬ 
canic  earthquakes  have  different  characteristics. 

Indicating  the  sum  of  the  frequencies  of  occurrence  over  an  amplitude 

range  from  infinity  to  some  value  A  by  N  ,  we  can  write: 

s 


N  =  /A  cA‘mdA 

S  'oo 


A1_in]^ 


(15) 

(16) 
(IT) 


In  practice  N  is  obtained  numerically  so  that  we  may  express;  it  in  summa- 

D 

tion  form: 


t  N;  — 

ff=» 


c 

1-m 


.1-m-iA 

A  L 


(13) 


The  term  b  in  Equation  .12  and  the  term  m  in  Equation  13  are  constants. 
The  physical  significance  of  these  constants  is  not  very  well  understood  and 
it  is  quite  remarkable  that  such  constants  exist.  In  the  investigation  of 
volcanic  quakes,  Minakami  has  found  that  m  has  different  values  depending  on 
the  type  of  quake  under  consideration. 
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A-type  volcanic  quakes  are  similar  to  tectonic  quakes  with  values  for 
ro  in  the  range  1.8  to  2.2.  E-type  volcanic  quakes  and  explosion  quakes 
seem  to  be  similar  and  have  much  larger  values  for  m  ranging  from  2.7  to 

4.0. 

Amplitude  data  for  the  Kilauea  Iki  eruption  were  analyzed  to  see  what 
sort  of  m  value  would  result.  When  the  data  are  plotted  as  log2.N  agaimt 
log  A,  the  curve  should  be  linear  and  have  slope  1-m  from  which  m  can  easily 
be  found.  Amplitude,  frequency  of  occurrence,  end  summation  values  are 
given  in  Table  2" .  The  data  are  shown  graphically  in  Figure  20  with  the 
value  of  m  equal  to  2.6. 

TABLE  20  AMPLITUDE  -  FREQUENCY  DATA,  KILAUEA  IKI 


Trace  Amplitude 

Frequency 

£n 

in  mm 

N 

1.0 

250 

342 

2.0 

42 

92 

3.0 

20 

50 

4.0 

6 

50 

5-0 

4 

24 

6.0 

4 

20 

7-0 

5 

l6 

8.0 

2 

n 

>8.0 

9 

9 

The  data  for  the  North  Pit  station  was  analyzed  separately  to  see  if 
this  method  might  reveal  something  unusual.  It  will  be  recalled  that  the 
North  Pit  data,  rarely  agreed  with  that  data  fiom  the  other  stations.  The 
value  of  m  was  found  to  be  2.0  which  can  be  considered  to  be  in  the  same 
range  as  the  value  determined  for  all  the  other  stations. 

A  similar  analysis  was  made  for  the  Vesuvian  quakes.  The  number  of 
data  points  was  smaller  and  the  curve  showed  a  much  different  character  - 
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two  distinct  slopes  are  apparent.  The  steeper  slope  gives  a  value  for 
m  =  4.8U,  which  would  place  Vesu1  us  in  the  explosive  type,  thus  ver¬ 
ily  ing  through  seismic  measurement  the  known  historical  character  of  Mount 
Vesuvius.  The  graph  of  this  data  is  shown  in  Figure  21. 

The  type  of  curve  resulting  from  the  graph  of  Vesuvian  amplitudes  has 
been  found  in  model  experiments  by  Mogi,  in  which  he  correlated  the  results 
with  the  mechanical  structure  of  the  medium,  which  in  turn  is  related  to 
tne  homogeneity  or  heterogeneity  of  the  medium. 

As  indicated  previously,  the  significance  of  the  coefficient  m  is  not 
toe  well  understood  at  this  time.  Perhaps  the  first  though*  that  arises  as 
a  possible  explanation  for  the  variation  in  m,  is  the  question  of  focal 
mechanism.  An  investigation  of  initial  motion  patterns  for  explosive  versus 
non-explosive  mechanisms  might  shed  some  light  on  this  problem.  Focal  depth 
may  be  a  significant  factor  since  the  m  values  for  A-type  and  B-type  quakes 
found  by  Minaivami  are  different.  Finally,  it  seems  quite  probable  that  the 
t;T?e  of  stress  application  may  be  siL  lificant.  In  the  case  of  a  tectonic 
quake,  the  stress  application  produces  a  mechanical  rupture  or  fault;  in  the 
case  of  a  volcanic  eruption,  the  stress  application  may  approximate  a  uni¬ 
form  internal  hydrostatic  system  giving  rise  eventually  to  an  explosive 
eruption.  In  the  case  of  an  externally  applied  intensive  stress  system, 
the  situation  is  again  different  and  may  be  revealed  by  a  characteristic 
value  for  m  just  as  the  tectonic  quakes  and  volcanic  explosion  quakes  have 
characteristic  values. 
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summary  ~AN~rr  cCttcLusioNs 


Seismic  effects,  +hat  is  earthquakes  occuring  in  volcanic  regions  and 
resulting  from  volcanic  processes,  have  been  the  subject  of  this  investi¬ 
gation.  Three  areas  selected  for  investigation  were  Italy,  Hawaii  and 
Japan.  Seismographs  from  these  areas  were  studied  and  statistical  analy¬ 
sis  was  applied  to  the  data.  The  following  main  conclusions  have  resulted 
from  the  investigation: 


1.  First  motion  patterns  can  usually  be  delineated  in 
non-explosive  volcanic  regions  which  indicate  that 
an  exclusive  first  motion  push  occurs  infrequently. 

2-  In  areas  where  the  volcanic  mechanism  is  explosive, 
a  first  motion  push  does  occur  but  not  exclusively 
and  many  other  seismic  events  also  occur  which  show 
a  quadrant al  pattern. 

3-  Interpretation  of  a  first  motion  push  is  subject 
to  serious  misinterpretation  if  location  of  the 
epicenter  and  focal  depth  are  unknown. 

4.  Volcanic  quakes  occur  in  swarms  and  are  generally 
of  low  energy  having  magnitudes  on  the  Richter 
Scale  that  range  from  plus  one  down  into  negative 
values . 
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5*  Volcanic  quakes  are  characterized  by  high  frequency 
motion  which  attenuates  rapidly  and  consequently 
are  not  usually  recorded  at  distances  beyond  a  few 
tens  of  kilometers . 

6.  Data  from  a  single  seismic  station  tends  to  show 
randomness  rather  than  reveal  identifying  char¬ 
acteristics.  Wen  used  in  conjunction  with  other 
stations,  the  data  then  responds  to  analysis. 

7*  The  frequency  of  occurrence  of  volcanic  quakes  shows 
promise  of  being  a  reliable  statistical  variable  with 

m  in  the  equation  NAm  =  c  being  the  important 
coefficient. 

8.  Tectonic  quakes  and  explosion  quakes  can  be  identi¬ 
fied  by  the  characteristic  m  value. 

9-  The  m  value  gives  promise  of  being  an  identifying 
characteristic  for  an  explosive  source  above  ground. 
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